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METHOD AND APPARATUS FOR COMPUTER MODELING OF AN ADAPTIVE 

IMMUNE RESPONSE 



COPYRIGHT NOTICE 

A portion of the disclosure of Ae patent document contains material that is subject 
to copyright protection. The copyright owner has no objection to the facsimile reproduction 
by anyone of the patent document of the patent disclosure, as it appears in the Patent and 
Trademark OfiBce patent file or records, but otherwise reserves all copyrigjit rights 
whatsoever. 

CROSS-REFERENCE TO RELATED APPLICATION 

The present invention is related to and claims priority xmder 35 U.S.C. §1 19(e) to 
U.S. Provisional Patent Application Serial No. 60/301,278, filed June 28, 2001, entitled 
**Method and Apparatus for Computer Modeling of T Cells," the specification of which is 
incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

The present invention relates generally to a computer model of the adaptive immime 
response. In one embodiment, the present invention relates to a computer model of the 
adaptive immune response within the firamework of signals conveyed at the site of antigen 
exposure, where the signals iuclude signals impacting antigen-presenting cells (APCs) and 
signals deUvered by APCs and by responding lymphocytes. In another embodiment, the 
present invention relates to an analytical model of an adaptive immune response. 

The human immune system has evolved as a complex process by which it is able to 
identify and respond to a vast array of genetically, biochemically, and behaviorally distiact 
microbial pathogens while not responding to the vast array of innocuous environmental 
elements or to the vast array of normal human cellular and biochemical elements. Adaptive 
immune responses involve populations of specialized immune cells or lymphocytes that 
have evolved to match the wide array of elements they may encounter. Particularly, 
lymphocyte populations are composed of a large number of individual cells, where each 
lymphocyte e^tpresses receptors with a distinct molecular sequence of a distinct afSnity and 
specificity. As a result, each lymphocyte only binds the molecular sequence, or antigen(s), 
that mblecularly interact with its receptor, and these cells are therefore classified as antigen- 
specific. 
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The nature of fhe adaptive immune response is to selectively activate and expand 
antigen-specific lymphocytes if and only if the antigens are presented in the correct context 
The dual requirements that an antigen-specific cell must not only recognize its particular 
antigen(s), but must also recognize the antigen within a particular enabling context, prevents 
5 lymphocyte responses to self or innocuous environmental antigens and promotes 

lymphocyte responses to pathogens. The enabling context is largely determined at the site of 
antigen exposure (i.e., peripheral tissue), which is generally anatomically separate from the 
site of lymphocyte expansion (i.e., secondary lymphoid tissue). 

APCs take up and process antigen at the site of antigen exposure and are primarily 

10 responsible for the transport of antigen from the site of antigen exposure to the site of 

lymphocyte expansion. Critical signals, directed to the APCs at the site of antigen exposure, 
shsqpe the manner in which APCs subsequently present antigen to lymphocytes. The manner 
of APC antigen pres^tation, in combination with other signals present in the secondary 
lymphoid tissue, is responsible for inducing a lymphocyte response and in guiding the 

IS character of the resfponse. For example, mediator production, costunulatory molecule 
expression, and antigen presentation by APCs varies according to the peripheral 
environment in which antigen is taken up and significantly contributes to the development 
of a T helper (Th) 1 or Th2 biased T lymphocyte response. 

Once activated, lymphocytes may directly or indirectly drive a cellular or humoral 

20 adaptive immune response, hi the case of an infecting pathogen, a cellular response is 

characterized by antigen-specific lymphocytes that traffic to the site of pathogen exposure. 
The lyn^)hocytes, which recognize antigen derived from the pathogen, may then directly 
kill the pathogen or activate other immune cells to kill the pathogen. A humoral response is 
characterized by antigen-specific lymphocytes that generate an antigen-specific antibody 

25 response; the antibodies bind the pathogens and facilitate their clearance from the body. As 
described above for Thl and Th2 lymphocyte responses, the generation of cellular and/or 
humoral adaptive immune responses is largely guided by a combination of antigen and 
context. 

The ad^tive immune response also generates memory lymphocytes, which allow 
30 the immxme system to increase its response efficiency. In the case of an infecting pathogen, 
these antigen-specific long-lived lymphocytes remain in the body after the pathogen is 
cleared, such that in later encounters with the same pathogen, the immune system responds 
more quickly and with greater strength than in the first encounter. Immune memory is 
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generated throughout the lifespan of an individual and confers the advantage that secondary 
exposure to a particular bacterium, virus, parasite, or fungus, can be cleared by the immune 
system with minimal compromise in the individual's ability to function. 

Because the unmime system must be able to respond to the vast array of pathog^ 
5 that may be encountered but should not respond to the even wider array of iimocuous 
environmental elements or self-elements, the adaptive immune response is a tightly 
regulated process. Hovrever, there is clearly a potential for inappropriate immune responses, 
as represented by the existence of allergic and autoiiomune diseases. 

The etiologies of inappropriate immime responses that manifest as allergic diseases 

10 (e.g., asthma, allergic rhinitis, food allergy) are improven but likely include genetic factors, 
history of exposure to environmental elements, and history of exposure to pathogens. The 
result is an inappropriate adaptive immune response to a normally irmocuous environmental 
• element (antigen), leading to elevated levels of inmumoglobulin (Ig) E and chronic 
inflammation at the exposure site. The development of phannaceutical treatments for these 

15 diseases has historically focused on controlling the symptoms of disease. However, as our 
understanding of immune processes improves, some newer treatments have been directed 
towards modifying the underlying inappropriate inumme response. This effort has been 
comphcated by the fact that the ad^tive response is highly complex, highly redundant, and 
tightly regulated, making the selection of appropriate intervention sites difficult 

20 The etiologies of inappropriate immune responses that manifest as autoimmune 

diseases (e.g., rheumatoid arthritis, multiple sclerosis, inflammatory bowel disease) are 
uiq)roven but likely include genetic factors, history of exposmre to environmental elements, 
and history of exposure to pathogens. The result is an inappropriate adaptive immune 
response to a self-molecule (antigen), leading to chronic inflammation and taigeted tissue 

25 destruction at the exposure site. As discussed above with allergic diseases, emerging 

phannaceutical therapies are directed towards modifying underlying inappropriate immune 
responses. However, the development of these new therapies has been complicated by the 
intricacies of the immune system and the need to selectively target inappropriate resfponses, 
while leaving appropxiztc immune responses intact. 

30 Several researchers have constructed simple mathematical models of 

antigen-specific lymphocyte expansion and its control by cytokines or antigen abundance 
(De Boer et al., J. ViroL, 75:10663-10669, 2001; Louzon ct al., /. Autoimmunity, 
17:311-321,2001; Yates etaL,y. TheorBiol, 206:539-200, 2000; Fishman&Perelson, 
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BuIL Math Biol, 61:403-436, 1999). These models were largely restricted to lyn^hocyte 
responses and did not represent important interactions that take place at the site of primary 
antigen exposure and that largely determine the enabling context for lymphocyte expansion. 
Specifically, these models did not include detailed representations of APC populations or 

5 the influence of particular peripheral tissues on the APCs. In addition, these models did not 
represent important feedback pathways from expanded lymphocyte populations to the site 
of antigen exposure; wherein, antigen-specific lymphocytes traffic from a lymphoid tissue 
to the site of antigen exposure to directly or indirectly act against the antigen source. 

Because these existing models do not include all aspects of the adaptive immune 

10 response, there is a need to develop a more comprehensive model of the adaptive immune 
response. One embodiment of the invention disclosed herein is a computer model of the 
ad^tive inomune response within the framework of signals conveyed at the site of antigen 
. exposure, where the signals include signals received by antigen-presenting cells (APCs) and 
signals delivered by APCs and by responding lymphocytes. 

15 

SUMMLARY OF THE INVENTION 

The present invention relates generally to a computer model of an adaptive immune 
response. One embodiment of the invention relates to a computer model of an adaptive 
immune response within the framework of signals conveyed at the site of antigen exposure. 

20 Another embodiment includes a representation of complex physiological regulatory 
mechanisms related to, for example, antigen-presenting cell (APC) recruitment, APC 
maturation, lymphocyte activation, and/or lymphocyte tra£E[cking. 

In one embodiment, the model can account for cellular dynamics and mediator 
production in response to antigen within a chronically inflamed peripheral tissue, as well as 

25 the regulatory effects on APC, APC population dynamics and activities, and lymphocyte 
population dynamics, including maturation, activation, effector fiinction, and apoptosis. hi 
addition, the model can account for inunune cell trafficking between a chronically inflamed 
peripheral tissue and secondary lymphoid tissues. In tiiis embodiment, the model can 
simulate a diverse set of adaptive immime responses, from acute to chronic progressive, and 

30 can predict the likely effects of therapeutic inventions. 

In another embodiment, a tolerant immune reaction can be modeled; wherein, 
signals and characteristics of the peripheral tissue do not provoke a lymphocyte response to 
a particular antigen. 
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Another CTibodiment of the invention is an analytical model of the adaptive immune 
response. 



BRIEF DESCRIPTION OF THE DRAWINGS 
5 FIG. 1 is a schematic representation of a computer system within which software for 

performing the methods of the invention may reside or be executed. 

FIG. 2 is a summary diagram that depicts the components of a computer model and 
their interconnectedness, according to an embodiment of the present invention where the 
airway is an example of a peripheral tissue. 
10 FIG. 3 depicts a flowchart for a method for developing a computer model of an 

adaptive immune response according to one embodiment of the invention. 

FIG, 4 depicts a flowchart for a method for developing a computer model of an 
adaptive immune response according to another embodiment of the invention. 

FIG. 5 A shows a portion of the Effect Diagram representing dendritic cell pC) 
IS precursors migrating from the blood into the airway and airway DC subpopulations 
undergoing maturation, according to an embodiment of the present invention. 

FIG. SB shows a portion of the Effect Diagram representing DC subpopulations 
trafBcking between the airway and lymph nodes and DCs in the lyn^h node, according to 
an embodiment of the present invention. 
20 FIG. 6A shows a portion of the Effect Diagram representing the biological processes 

involved in migration and maturation of airway DCs, according to an embodiment of the 
present invention. 

FIG. 6B shows a portion of the Effect Diagram representing the biological processes 
involved in maturation and migration of lymph node DCs, according to an embodiment of 
25 the present invention. 

FIG. 7A shows a portion of the Effect Diagram representing the biological processes 
involved in antigen presentation by airway DCs, according to an ^bodiment of the present 
invention. 

FIG. 7B shows a portion of the Effect Diagram representing the biological processes 
30 involved in antigen presentation by lymph node DCs, according to an embodiment of the 
present invention. 
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FIG. 8 shows a portion of the Effect Diagram representing the biological processes 
of mediator production by airway and lymph node DCs, according to an embodiment of the 
present invention. 

FIG. 9A shows a portion of the Effect Diagram representing the biological processes 
5 that determine IL-12 homodimer and heterodimer production capability of maturing airway 
DCs, according to an embodiment of the present invention. 

FIG. 9B shows a portion of the Effect Diagram representing the biological processes 
that determine IL-12 homodimer and heterodimer production capability of DCs trafficking 
between the airway and lymph nodes and maturing lymph node DCs, according to an 
10 embodiment of the present invention. 

FIG. 9C , in conjunction with FIG. 9D, shows a portion of the Effect Diagram 
representing-the biological processes involved in late-stage regulation of IL-12 homodimer 
and heterodimer production capability in mature DCs, according to an embodiment of the 
preset invention. 

15 FIG. 9D, in conjunction with FIG. 9C, shows a portion of the EflFect Diagram 

. . representing tiie biological processes involved in late-stage regulation of IL-12 homodimer 
and heterodimer production capability in mature DCs, according to an embodiment of the 
. present invention. 

FIG. 10 shows an example of an Effect Diagram that displays a taxonomy of CD4+ 
20 T lymphocyte states, according to an embodiment of the present invention. 

. FIG. 11 shows an example of a portion of an Effect Diagram that calculates CD4+T 
lyn^hocyte cytokine production, according to an embodiment of the present invention. 

FIG. • 12 shows an example of a portion of an Effect Diagram relating cytokine 
binding and cell-cell interactions to CD4+ T lymphocyte exp^ion, differentiation^ and . 
25 apoptosis, according to an embodiment of the present invention. 
I FIG. 13 shows the Effect Diagram of FIG. 12 with an example of the mathematical 

calculation contained in the LN pe2 expansion function node, according to an embodimdnt 
of the present invention. 

FIG. 14 illustrates an Effect Diagram that depicts Thl and Th2 cell trafficking as 
30 modulated by chemokines and adhesion molecxiles according to an embodiirient of the 
present invention where the airway (AW) is an example of a peripheral tissue. 
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FIG. 15 illustrates an output of the model showing expansion of Thl (♦) and Th2 
cells (•) as a function of antigen dose after simulated in vitro primary culture, according to 
an embodiment of the present invention. 

FIG. 16 illustrates an output of the model fhat depicts evolving numbers of Thl (-) 
5 and Th2 (— ) cells during a simulated in vitro primary culture at relatively low antigen dose, 
according to an embodiment of the present invention. 

FIG. 17 illustrates an output of the model that depicts evolving numbers of Thl (-) 
and Th2 (— ) cells during a simulated in vitro primary culture at a lower antigen dose than 
used in FIG. 1 6, according to an embodiment of the present invention. 
10 FIG. 1 8 illustrates an Effect Diagram that depicts the recruitment of DC precursor 

populations from flie blood into the peripheral tissue according to an embodiment of the 
present invention where the airway is an exanq)le of peripheral tissue. 

FIG. 19 illustrates an output of the model that depicts the simulated kinetics of 
adhesion molecule CTcpression in the airway tissue following antigen challenge at time equal 
IS to ^zero, according to an embodiment of the present invention. 

FIG. 20 illustrates an output of the model that depicts the simulated kinetics of blood 
monocytes (— )compaied with the experimental results reported by Whitelaw 1966 (•), 
according to an embodiment of the present mvention. 

FIG. 21 illustrates an output of the model that depicts the simulated kinetics of tissue > 
20 DCs (— ) compared with the experimental results reported by Holt et al. 1994 (•), according 
to an embodiment of the present invention. 

FIG. 22 illustrates an output of the model that depicts the kinetics of blood 
monocytes (-), blood DC (- and lung DC (- ■ — ) responses following antigen 
challenge at time equal to zero compared with the experimental results reported by Upham 
25 et al. 1999 (blood monocytes = o; blood DCs = x) and McWilUam et al. 1994 (lung DCs = 
+), according to an embodiment of the present invention. * 

FIG. 23 illustrates an output of the model that depicts the sunulated kinetics of DC ' 
migration to the lymph node in the context of non-productive interactions with CD4+.T 
lymphocytes (— ) compared with the experimental results reported by InguUi et al. 1997 (•) • 
30 and Veimaelen et al. 2001 (■), according to an embodiment of the present invention. 

FIG. 24 illustrates an output of the model that depicts the simulated kinetics of DC 
migration to the lymph node in the context of productive interactions with CD4+ T 
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lymphocytes compared wifli the experimental results reported by Ingulli et al. 1997 (•), 
* according to an embodiment of the present invention. 



DETAILED DESCRIPTION 

Overview 

The present invention relates to computer modeling of an adaptive immune 
response. The adaptive immune response model can be used in isolation or integrated with 
other components to represent a healthy or diseased physiological system whose state is 
afiFected by the adaptive immune response. Embodiments of the present invention relate to 
modeled responses of antigen-presenting cells (APCs) and lymphocytes to immunogenic 
stimuU in the context of human diseases that involve the adaptive inamvme response (e.g., 
allergic asthma). In particular, one embodiment of the model includes a peripheral tissue 
environment, a lymphoid tissue environment, and traffic of immune cells between the two 
compartments. The peripheral tissue represented in the model can include for example, 
lung, skin, intestine, joint, or the central nervous system. The term lymphoid tissue 
environment as used herein can include primary, secondary, and tertiary lyn5)hoid tissues.* 

The model can further include the character and kinetics of antigen exposure, the 
character and dynamics of APC populations (e.g., dendritic cells), ttie expansion, 
. differentiation and pontraction of antigen-specific lymphocyte populations (e.g., T 
lymphocytes), and the creation and maintenance of memory lymphocyte populations. 

In one embodiment, the model includes one or more of the following features: (1) 
comrnunication between a peripheral tissue and a secondary lymphoid tissue, (2) 
establishing a< stable balance of popiilations within and between the two compartments, and 
(3) feedbadc.pathways that allow for progressive changes in a disease state (i.e., the model 
is not limited to a steady-state representation). 

The conq)uter model of the present invention can be used to identify pharmaceutical 
•interventions to treat immune diseases such as allergic asthma. In another embodiment, 
therapies affecting pathways that are present in die model can be implemented and used to 
•predict therapeutic outcomes. 

Mathematical Model 
The mathematical model implemented by the computer-executable software code 
represents the dynamic biological processes related to an adaptive immune response. The 
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fonii of the mathematical equations employed may include, for example partial differential 
equations, stochastic differential equations, differential algebraic equations, difference 
equations, cellular automata, coupled maps, equations of networks of Boolean or fuzzy 
logical networks, etc. In one embodiment, the forms of the mathematical equations used in 
5 the model are ordinary differential equations : 

dx/dt = f(x,p,t), 

where x is an N dimensional vector whose elements represent the biological variables of the 
system (for example concentrations of chemical mediators, monocyte density in the blood, 
naive T cell density in lymphoid tissue, etc.), t is time, dx/dt is the rate of change of x, p is 
10 an M dimensional set of system parameters (for example sensitivity of the blood monocyte . 
recruitment rate to P-selectin, naive T cell efflux rate, equiUbrium dissociation constant for 
IL-IQ, etc; ), and f is a function that represents flie complex interactions among biological - 
variables. • . 

The term '^biological variables" refers to the biological constitueiits that make up a 
IS : biological process. Mathematically, the above x rq>resents the biological variables in the 
. model. For example, the biological variables can include metabolites, DNA, RNA, proteins, 
enzymes, hormones, cells, organs, tissues, portions of cells, tissues, or organs, subcellular 
organelles, chemically reactive molecules like H^, superoxides, ATP, citric acid, protein 
albumin, as well as combinations or aggregate representations of these types of biological 
20 variables. In addition, biological variables can include response-provoking agents, such as • 
antigen or mefhacholine, and tfaer^eutic agents such as steroids, P-agonists, or leukolriene 
antagonists. 

The term ''parameter'' is used herein to mean a number that characterizes the 
behavior of a single biological variable or the interaction between two or more biological 

25 variables. For example, a parameter could be the baseline synthesis of a mediator, baseline 
expression of a cell surface molecule, or the maximum number of lymphocytes that may 
interact with any one APC. Parameters may also be used to specify synthetic or 
environmbntal factors, as well as intrinsic biological properties. 

The term **biological process" is defined herein to mean an mteraction or series of 

30 interactions between biological variables. Biological processes can include, for example, 

cellular recruitment, regulation of maturation, induction of anergy, or regulation of chemical 
mediator production. Each biological variable of the biological process can be influenced. 



9 



wo 03/001891 PCT/US02/20672 

for example, by at least one other biological variable in the biological process by some 
biological mechanism, which need not be specified or even understood. 

The term *T>iological state" is used herein to mean the result of the occurrence of a 
series of biological processes. As the biological processes change relative to each other, the 
5 biological state also imdergoes changes. One measurement of a biological state, is the level 
of activity of biologic variables, parameters, and/or processes at a specified time and under 
specified experimental or environmental conditions. 

• In on6 embodiment the biological state can be mathematically defined by the values 
of X and p at a given time. Once a biological state of the model is mathematically specified, 
10 numerical integration of the above equation using a computer determines, for example, the 
time evolution of the biological variables x(t) and hence the evolution of the biological state 
• overtime. 

A biological state can include, for example, the state of an individual cell, a 
population of cells, a tissue, and/or a multi-cellular organism: A biological state can also 
15 include the state of a mediator concentration in the plasma, interstitial fluid, or intracellular 
fluid. :For example, a biological state ofanAPC population can include the APC density or 
the antigen-presenting capacity in a particular peripheral tissue of a particular patient type at 
. a particular point in time. 

Ad£q)tive immime responses are the set of responses that target specific antigra for 
20 ; immune activity and that are mounted by the immune system. In one embodiment of the 
invention, the biological state modeled is the state of an adaptive immune response. The 
tenn "adaptive immune response" as used herein comprises a combination of at least two of 
.the following classes: biological processes at the site of antigen exposure, the impact of 
these biological processes on the character or behavior of immune cells, biological 
25 processes relating to cellular dynamics between a site of antigen exposure and lymphoid 
. tissue, biological processes relating to a primary lymphoid tissue, biological processes 
relating to the interaction between APCs and lymphocytes, or biological processes relating 
to the feedback of immune cells on biological processes at the site of antigen exposure. 

For example, the site of antigen exposure can include a peripheral tissue where the 
30 . biological processes of the tissue can include invasion of host cells by a pathogen, 
interactions among inflammatory cells, or mediator levels. The impact of biological 
processes on immune cells might include, for example, the impact of gamma interferon 
production and binding on the ability of APCs to produce interleukin 12. Cellular dynamics 

10 
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between a site of antigen exposure and a lymphoid tissue might include for example, 
trafficking of APCs through the lymphatics or trajfficking of lymphocytes through the blood 
Feedback of immune cells on biological processes at the site of antigen exposure might 
include, for example, T lymphocyte activation of peripheral tissue macrophages. In one 
embodiment, a model of the adaptive immune response includes representation of at least 
two biological compartments including a peripheral site of antigen exposure and a lymphoid 
tissue, representation of the biological processes at a peripheral site and their impact on 
APCs and lymphocytes, representation of the cellular dynantiics between a peripheral site 
and a lymphoid tissue, and representation of the impact of immune cells on a peripheral site. 

The model of the adaptive immune response could be integrated with any number of 
other componmts to represent a healthy or diseased biological state. The method of this 
integration would involve modifying fho peripheral tissue to assume the attributes of the 
peripheral tissue targeted in a particular biological state. In one embodimmt, characteristics 
of cell types, cellular abnormalities, physiological abnormalities, and tiie diemical mediator 
environment of the peripheral tissue can be modeled qipioprialely. In another embodiment, 
the adaptive imoiune response model can be modified to reflect the nature of APCs and 
lymphocytes tiiat associated with a biological state. The regions of inter&ce can, for 
example, include modulation of antigen-presenting cell function by the peripheral tissue and 
modulation of the peripheral tissue by the antigen-presenting cells and the antigra-specific 
lymphocytes. 

The tenn "sunulation" is used herein to mean solution of a mathematical model by 
the numerical or analytical methods. For example, simulation can mean the numerical 
integration of the mathematical model of the biological state defined by the above equation, 
(i.e. dx/dt = f(x, p, t)) and specifying an initial value of x. 

The term "disease state" is used herein to mean a biological state where one or more 
biological processes are related to the cause(s) or the cUnical signs of the disease. For 
example, a disease state can be the state of a diseased cell, a diseased organ, a diseased 
tissue, or a diseased multi-cellular organism/ Such diseases can include, for example, 
acquired immune deficiency syndrome, delayed-type hypersensitivity, systemic 
anaphylaxis, allergic asthma, cancer, inflammatory bowel disease, systemic hxpvis 
erythematosus, multiple sclerosis, type I diabetes, and rheumatoid arthritis. A diseased 
multi-cellular organism can be, for exanq>le, an individual patient, a specific groiq> of 
human patients, or flie genenil human population as a whole. A diseased state could refer to. 
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for example, a diseased protein such as a defective interferon-gamma receptor or a diseased 
process, such as defects in cellular activation, cell signaling, or cell mediator production, 
which may occur in several different organs. 

The term **biological attribute" is used herein to mean observations or diagnostic 
criteria associated with a biological state. The biological attributes of a biological state can 
be measurements of biological variables, parameters, and/or processes. For example, for 
the disease allergic asthma, the biological attributes associated with the adaptive immune 
respoiise can include APC dynamics, T lymphocyte dynamics, or T lymphocyte cytokine 
production. . 

The term "simulated biological attribute" is used herein to mean measurements on 
-model variables or processes corresponding to biological attributes. For example, simulated 
hiological attributes associated with tiie modeled adaptive iminune response might -include 
measurements of APC or T lymphocyte dynamics. 

: The teim ^^substantially consistenf * is used herein to mean fliat the relationship 
between simulated biological attributes and biological attributes is sufficiently similar to 
;conclude that a simulated biological attribute accurately represents a biological attribute; tiie 
sunulated biological attribute and biological attribute do not have to be identical. The term 
'"substantially consistenf ' can be, for example, simulation outcomes demonstrating relative 
changes in a pattem of cytokine levels tiiat are similar to relative changes in a pattern of 
cytokine levels measured in an in vitro e3q)eriment but with different absolute values. 

The term "Reference pattem" is used herein to mean a set of biological attributes that 
are measured in a normal or diseased biological system under spedjSed experimental 
conditions. For example, the reference pattem of an aUergic asthmatic might include 
measurements performed on lung exudate via broncho-alveolar lavage, lung function via 
spirometry, lung tissue via biopsy, or blood via venipuncture at a specified time following a 
. particular chemical mediator or antigen stimulus. Alternatively, the reference pattem of 
. APC behavior might include measurements on cell cultures derived fiom a normal or 
diseased human or animal under defined conditions. 

Computer Svstem 

FIG. 1 shows a system block diagram of a computer system within which the 
methods described above can operate via software code, according to an embodiment of the 
present invention. The computer system 1(K) includes a prcx:essor 102, a n^ 
and a static moDOLory 104, which are coupled by bus 106. The computer system 100 can • 
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further include a video display unit 108 (e.g., a liquid crystal display (LCD) or cathode ray 
tube (CRT)) on which a xiser interface can be displayed The compute: system 100 can also 
include an alpha-numeric input device 110 (e.g., a keyboard), a cursor control device 1 12 
(e.g., a mouse), a disk drive unit 1 14, a signal generation device 116 (e.g., a speaker) and a 
network interface device medium 118. The disk drive imit 114 includes a computer- 
readable medium 1 1 S on which software 120 can be stored. The software can also reside, 
completely or partially, within the main memory 103 and/or within the processor 102. The 
software 120 can also be transmitted or received via the network interface device 118. 

The temi "computer-readable medium" is used herein to include any medium which 
is capable of storing or encoding a sequence of instructions or codes for performing the 
methods described herein and can include, but not limited to, optical and/or magnetic 
storage devices and/or disks, and carrier wave signals. 

The Compater Model 
• •* Suitably, a con^^uter model can be used to inclement at least some embodimmtso^^ 
the presmt invention. The computer model can be used for a variety of purposes. For 
sample, it can enable a researcher to: (1) simulate flie dynamics of ttie biological state 
associated with an adaptive immune response, (2) visualize key biological pathways for the 
initiation and maintenance of an adaptive immune response and the feedback within and 
between these pafliways, (3) gain a better understanding of the physiology of an ad^^tive 
immune response, (4) explore and test hypotheses about adaptive immune responses, (5) 
identify and prioritize potential therapeutic targets, (6) identify dififerent types of response 
and their underlying mechanisms, (7) identify surrogate maikors of response types, and (8) 
organize knowledge and data that relate to the adaptive immune response. 

In addition to simulation c^abUities, the computer model can include a built-in 
database of references to the scientific literature on which the model is based. Users can 
augment this database with additional references or other comm^tary and can link the 
informatibn to the relevant component. The computer model can be a multi-user systrai in 
which the information can be shared throughout an organization. Thus, the computer model 
can be a specialized knowledge management system focused on the adaptive immune 
response. . 

While fbc following discussion is m terms of a computer model, one of skill in the 
art would recognize that the mathematical equations of the model may be analytically or 
numerically inqilemented widiout the assistance of a compvtier. 
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Effect Diagram and Summary Diagram 

In one embodiment, the computer model can represent various biological 
components or mechanisms through the use of an Effect Diagram, including a smnmary 
diagram and more detailed modules that represent the various biological processes of the 
biological system being modeled. These modules provide not only a conceptual map of the 
model, but also represent and encode sets of ordinary differential equations for numerical 
integration, as discussed more fully below in the section entitled 'Mathematical Equations 
l^coded in the Effect Diagram**. ; 

FIG. 2 is a summary diagram that depicts the adsqptive immune response components 
' of a model and flieir interconnectedness, according to an embodiment of the present 
invention, where the airway is an example of a peripheral tissue. Each squared node shown 
in FIG. 2 represents a functional module diagram discussed more fully below. The 
combined functional Effect Diagrams can represent and model flie lecniitment, phenotyplc 
jnatuiation, antigen processing, death, and departure of APCs in a iq>resentative peripheral 
tissue and lymph node (LN); the recniitment, expansion, differentiation, death, and 
departure of T lymphocytes in a representative LN and pmpheral tissue; and the produdtion 
of chemical mediators by these cells, subject to tibie time-varying stimuli of antigen and 
inflammatory signals fiom the peripheral tissue. 

In one embodiment, antigen presentation can be characterized in the computer 
model by the availability of APCs, specifically dendritic cells (DCs) over time. Antigm 
presentation is fiirther characterized by the average antigen density per cell, the average 
expression of costimulatory molecules (e.g., CD80, CD86) per cell, and DC mediator 
production . The degree of costimulation accompanying antigen presentation is also 
determined by the average expression of costimulatory counter-receptors (e.g., CD28) on T 
lynophocytes being activated. Costimulation effects involving other accessory molecules- 
can be modeled indirectly. For example, CTLA-4 can be modeled indirectly by varying the 
effectiye role of CD28. 

T lymphocyte population dynamics in the mathematical model are regulated througih 
cell-cell interaction, cytokine production, and cytokine effects. In each representative tissue 
compartment, cytokme producing cells contribute to a conunon cytokine * 
. subpopiiMonsofT lymphocytes -described below -reqK)^ 
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the combined T lymphocyte populations and the ambient cytokine milieu vary over time in 
tand^, each influencing the other. 

The particular Effect Diagrams shown below are discussed in reference to particular 
biological functions (e.gi, DC recruitment; DC states, T lymphocyte states). Pages A-1 
through A- 17 show the complete set of Effect Diagrams included in the present embodiment 
of the invention. • 



Mathematical Equations Encoded in the Effect Diagram 

As mentioned above, the EfTect Diagram is a visual representation of the model 
equations. This section describes how the diagram encodes a set of ordinary differential 
equations. Note that altihough die discussion below regarding state and function nodes 
refers to biological variables for consistency, tfie discussion also relates to variables of any 
appropriate type and need not be linuted to just biological variables. 
State andFunction Nodes 

State and function nodes display the names of ttie variables ttiey represent and their 
location in the model. Their arrows and modifiers indicate their relation to oUbet nodes 
within the model. State and function nodes also contain the parameters and equations that 
are used to compute the values or their variables in simulated experiments. In one - 
embodiment of the computer model, the state and function nodes are generated according to 
the method described in US. Patrat 6,051,029 and co-pending q)plication 09/588,855, both 
of which are entitled 'Method of generating a display for a dynamic simulation model 
utilizing node and link represoitations," and both of which are incorporated herein by 
reference. Further exainples of state and fimction nodes are further discussed bel^^^ 

State nodes, the single-border ovals in ttie Effect Diagram, represent 

State Node ) 

' variables in the system the values of which are determined by the cumulative 



effects of Its inputs over time. 

State node values are delBned by differential equations. The predefined parameters 
for a state node include its initial value (So) and its status. State nodes that have a half-life 
have the additional parameter of a half-life (h) and are labeled with a half-life i^symbol. 

Function nodes, tiie double-border ovals m the Effect Diagram, 
represent variables in the system the values of wfaidi, at any point in time, are 
detemuned by inputs at that same point in time. 
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Function nodes are defined by algebraic functions of their inputs. The predefined 
parameters for a fiinction node include its initial value (Pq) and its status. 

Setting the status of a node effects how the value of the node is determined. The status 
of a state or function node can be 

• Computed — the value is calculated as a result of its inputs 

• Specified-Locked — the value is held constant over time 

• Specified Data — the value varies with time according to predefined data points. 

State and function nodes can appear more dian once in the Bffect Diagram as alias 
liodes. Alias nodes are indicated by one or more dots, as in the state node illustration above. 
All nodes are also defined by their position, with respect to arrows and other nodes, as being 
either source nodes (S) or target nodes (T). Source nodes are located at the tails of arrows, 
and target nodes are located at the heads of arrows. Nodes can be active or inactive. Active 
nodes are white. Inactive nodes match the backgroimd color of the Effect Diagram. 
State Node Equations 

The computational status of a state node can be Computed, Spedfied-Locked, or 
Specified Data. 



State Node Computed ^ =y 

dt 



sumof arrov^terms y/Aen h^O 

—^S(f) + sumqf arrowierms when A > 0 



Where iS is the node value, t is time, S(t) is the node value at time, t, and h is the half-life. . 
The three dots at the end of the equation indicate there are additional terms in the equation 
resulting fix>m any effect arrows leading into it and by any conversion arrows that lead out 
of it. If h is equal to 0, then the half-life calculation is not performed and dS/dt is determined 
solely by the arrows attached to the node. 

State Node Spscifiedr Locked S(t) ^S^jor^aUi 

State Node Specified Data S(t) is defined specified data entered for flie state node. 

State node values can be limited to a minimum value of zero and a maximum value 
of one. If lindted at 2m>, iS can never be less than zero and the valiie for 5 is m 
it goes negative. If limited at one, S cannot b e greater flian one and is reset to one if it 
exceeds one. 
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Function Node Equations 

Function node equations are computed by evaluating the specified function of the 
values of the nodes with airows pointing into the function node (arguments), plus any object 
and Effect Diagram parameters used in the function expression. To view the specified 
function, cUck the Evaluation tab in the function node Object window. 
The Effect Diagram — Arrows 

Arrows link source nodes to target nodes and represent the mathematical relationship 
between the nodes. Arrows can be labeled with circles that indicate the activity of die auow. 
A key to the annotations in the circles is located in the upper left comer of each module in 
the Effect Diagram. If an arrowhead is solid, the effect is positive. If the arrowhead is 
hollow, the effect is negative. 

Arrow Types 

._^^^^^y^Sect arrows, the thin arrows on the Effect Diagram, link source state or function 

f nodes to target state nodes. Effect arrows cause changes to target nodes but have no 
effect on source nodes. They are labeled with circles ttiat indicate the activity of the arrow. 

-O 

Conversion arrows, the thick arrows on the Effect Diagrami, represent the way the 
. contents of state nodes are converted into the contents of the attached state nodes. They are 
labeled with circles that indicate the activity of the arrow. The activity may effect fhe source 
node or the target node or both nodes. The conversion can go either way. 

4 1^ Argument arrows specify which nodes are input arguments for function nodes. They 

do not contain parameters or equations and are not labeled with activity circles. 

Arrow Characteristics 

Effect or conversion airows can be constant, proportional, or interactive. 

— li^Arrows that are constant have a break in the arrow shaft. They are used when thfe rate 
25 of change of the target is independent of the values of the source and target nodes. 

Anows fliat are proportional have solid, unbroken shafts and are used when the rate 
of changcis d^mdent on, or is a function oi^ fhe v^hies of the source node. 
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^Arrows that are interactive have a loop from the activity circle to the target node. 
.-i/A|to.They indicate that the rate of change of the target is dependent on, or a function of, 
the value of bo& the source node and the target node. 

Arrow Properties can be displayed in an Object window (not shown). The window may 
5 also include tabs for displaying Notes and Arguments associated with the auow. If Notes 
are available in the Object window, the arrow is labeled with a red dot (•). 

Arrow Equations: EflFect Arrows 

"Proportional Effect Arrow: The rate of change of target tracks source node value. * 
dt 

10 Where 7is the target node, C is a coefficient, S is the source node, and a is 

an exponent. 

Constant Effect Anew: The late of change of the target is constant 

dt 

Where T is the target node and is a constant. 

IS Interaction Effect Arrow: The rate of change of the target depends on both the source node 
and target node values. 

€11 

Where T is the target node, S is the source node, and a and b are exponents. 
' This equation can vary depending on the operation selected in the Object 

20 window; The operations available are S+r, S-Z S*T, T/S, and S/T. • 

Arrow Equations: Conversion Arrows 

Proportional Conversion Arrow: The rate of change of the target tracks flie value of source 
node. 

dT 

25 Where T is the target node, S is the source node, C is a coefficient, /{ is a 

conversion ratio, and a is an exponent. 

Constant Conversion Arrow: The rates of change of target and source are constant such 
that an increase in target con:espQiids to a decrease in source. 
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dt 

d£ 

Where J is the target node, S is the source node, iT is a constant, and i2 is a 

conversion ratio. 

Interaction Conversion Arrow: The rates of change of the target and source depend on both 
source and target node values such that an increase in target corresponds to a decrease in 
source. 

Where T is the target node, iS* is the source node, a and b are exponents, and 
is a conversion ratio. This equation can vary dqiending on the operation 
selected in fhe Object window. The operations available are iS+ T, S-Z S*T, 
m, and iS/T. 

The Effect Diagram — Modifiers 

Modifiers indicate the efifects nodes have on the arrows to which they are connected. The 
type of modification is qualitatively indicated by a symbol in the box. For example, a node 
can allow block [U, regulate (El, inhibit E], or stimulate [Elan arrow rate. 

A key to the modifier annotations is located in the upper left comer of each module. 

Modifier Properties can be displayed in the Object Window. The window may also include 
tabs . for displaying tiie notes, arguments, and specified data associated with the modifier. If 
notes are available in the Object window, the modifier is labeled with a red dot («). ( 



dT Ju\ 



Effect Arrow, Modifier Equation : 
Where T is the target node, M is a multiplier constant, iV is a normalization constant, f()is^ 
function (either linear or specified by a transform curve), and arrowterm is an equation 
fi^gment fiiom the attached arrow. 

Modifier ftffect 
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By default, conversion arrow modifiers affect both the source and target arrow teims. 
However, in some cases, a unilateral, modifier is used. Such modifier will ajBfect eitha: a 
source arrow term or on target arrow term; it does not affect both arrow terms. 

Conversion arrow, Source Only Modifier Equation: 

• — = M • / — • arrowterm + other attachedarrow terms 

Conversion arrow, Target Only Modifier Equation: 
— = M • /I — I •arro\9tenn + other attached arrow terms 

The equation for a source and target modifier uses both the Source Only equation and the 
Target Only equation. 

When mxiltiplicative and additive modifiers are combined, effect is given precedence. For 
example, if the following modifiers are on an arrow, 

al,a2: Additive, Source and Target 
ml,m2: Multiplicative, Source and Target 

A1,A2: Additive, Target Only 

Ml JM2: Multiplicative, Target Only 

then the rates 'are modified by 

Target node: (al+a2+Al+A2) » (ml*m2) * (M1*M2) 
Source node: (al+a2) * (ml''ii2) 



Emhnriiiii ents of the Invention 

no. 3 depicts a flowchart for a method for developing a computer model of an 
ad2Q)tive immune response according to one embodiment of the invention. At step 3 10, data 
relating to a biological state oftiiead^tive immune response is identified At step 320, 
biological processes related to the data are identified. These biological processes define at 
least one portion of tiie biological state of the adaptive immune response. At step 330, the 
biological processes are combined to form a simulation of the biological state of the 
adsqptive immune response. 

Another embodiment of the invention is a method of developing a computer model 
of the adaptive immune response wherein the biological state of the adaptive immune 
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response is a biological state of an acute response. In another embodiment of the invention 
is a method of developing a computer model of the adaptive immune response wherein, the 
biological state of ttie adaptive immune response is a biological state of a chronic response* 
In another embodiment, at least one biological process is associated with a biological 
variable that is a therapeutic agent, A therapeutic agent of the invention can be, for 
example, recombinant IL-12, TNF-alpha blockade, steroids, or phosphodiesterase inhibitors. 

The method for developing a computer model of an adaptive immime response can 
further coiiq)rise the optional steps of 340, 350, 360, and 370 for validating the computer 
model, as depicted in FIG. 3 . In the validation process, at step 340 a simulated biological 
attribute associated with the biological state of the adaptive immune response is produced. 
At step 350, the simulated biological attribute is compared with a correspondiug biological 
attribute in a reference pattern of the adaptive immune response. At steps 360 and 370, the 
validity of the conq>uter model is identified. At step 360, it is detemiined whether the 
simulated biological attribute is substantially consistent with the biological attribute 
associated with the reference pattern of the adaptive immune response. At step 370, if the 
simulated biological attribute is substantially consistent with the biological attribute 
associated with the reference pattern of the ad^tive immune response the computer model 
is identified as a valid computer model of an adaptive immune response. 

FIG. 4 depicts a flowchart for a method for developing a computer model of an 
adaptive immune response according to another embodiment of the invention. At step 410, 
data relating to a biological state of the adaptive immune response is identified. At step 
420, biological processes related to the data are identified. These biological processes 
define at least one portion of the biological state of the adq}tive immune response. At step 
430, a first mathematical relation among biological variables associated with a first 
biological process fiom the biological processes is fonned. At step 440, a second 
mathematical relation among biolo^cal variables associated with the first biological process 
and a second biological process associated with the biological processes is fonned. 

Steps 450, 460, and 470 can be optionally performed to produce a simulated 
biological attribute that is substantially consistoit with at least one biologplcal attribute 
associated with a reference pattern of the adaptive immune response. At conditional step 
450, a detennination is made as to whether a simulated biological attribute or a series of 
simulated biological attributes is to be produced. If a simulated biological attribute is to be 
produced, &e process continues to step 460. At step 460, a set of parametric changes in the 
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first mathematical relation and the second mathematical relation is created. At step 470, a 
simulated biological attribute based on at least one parametric change from the set of 
parametric changes is produced. 

Steps 480, 490, SOO, SIC, and S20 can be optionally performed to obtain a 

5 representation of the chronological progression of a diseased adaptive immune response, for 
example from, a healthy state to a disease state. At step 480, a determination is made as to 
whether a biological variable or a parameter is converted. If a biological variable is to be 
converted the process proceeds to steps 510, and 520. At step 5 1 0, a first biological 
variable is converted into a converted biological variable the value of which changes over 

10 jtime. This first biological variable is associated with at least one from the first 

mathematical relation and the second mathematical relation formed in steps 430 and 440. 
At step 520, a series of simulated biological attributes are produced based on the converted 
biological variable. The series of simulated biological attributes are substantially consistent 
with a corresponding biological attribute associated with a reference pattern of the adaptive 

15, immune response. The series of simulated biological attributes represent the chronological 
progression of corresponding biological attributes in the reference pattern of the adaptive 
immune response. If a parameter is to be converted to obtain a series of simulated 
' biological attributes, the process proceeds to steps 490 and 500. At step 490, a parameter is 
converted into a new biological variable the value of which changes over time. This 

20 parameter is associated with at least one from the first mathematical relation and the second 
mathematical relation formed in steps 430 and 440. At step 500, a series of simulated 
biological attributes are produced based on the converted biological variable. 

Another embodiment of the present invention is a method for developing a computer 
model of an adaptive immune response fliat includes the steps of identifying data related to 

25 the biological state of the adaptive immune response; identifying biological processes 

related to the data, the biological processes defining at least one portion of the biological 
* state of the adaptive immune response; and combining the biological processes to form a 
simulation of the biological state of the adaptive immune response in tibie context of a 
per^heral tissue environment and a lynq)hoid tissue enviromnent 

30 Another embodiment of ttie invwtion is a method of developing a computer model 

.of the adq[>tive immune response wherein at least one biological process is associated with 
•xecniitmeDt of immune cells into ttie peripheral tissue environment In a fiirOier 
embodiment, the method includes immune cells that are blood dendritic cells and blood 
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monocytes. In yet a further embodiment, the method of developing a computer model of an 
adaptive immune response includes the combining of biological processes so that the 
peripheral tissue environment is modeled with preferential recruitment of the blood 
dendritic cells over the blood monocytes. 

Another embodiment of the invention is a computer model of a biological state of an 
adaptive immune response. The computer model comprises Code to define biological . 
processes related to the biological state of the adaptive immune response; and code to define 
mathematical relationships related to interactions among biological variables associated with 
the biological processes. At least two biological processes fiom the biological processes are 
associated with the mathematical relationships. The combination of the code to define the 
biological processes and the code to define mathematical relationships define a simulation of 
the biological state of flie ad2q)tive immune response in the context of a peripheral tissue 
environment and a lymphoid tissue environment. 

. In one embodiment oftfaeconQ>ut^ model ofthead^ttve immune response at le^ 
one biological process is associated with recruitment of immune cells into the peripheral 
tissue environment In a fiirther embodiment, the model includes immune cells that are blood 
dendritic cells and blood monocytes. In yet a fiutfaer embodiment, the computer model of an 
adaptive immune response includes die combining of biological processes so that the 
peripheral tissue environment is modeled with preferential recruitment of the blood dendritic 
cells over the blood monocytes. 

The computer model can further comprise code to define two compartm^ts, wherein 
one compartmmt includes biological processes related to a peripheral tissue ravironment 
and the second compartment includes biological processes related to a lymphoid tissue 
environment Further, die computer model can include a code to define the interaction 
between these two compartments. 

Yet another embodiment of the invention is a computer executable software code that 
comprises of code to define biological processes related to a biological state of an adaptive 
immune response including code to define mathematical relations associated with the > 
biological processes. The biological processes defined by the code are associated with the 
biological state of the adaptive immune response. 

The computer executable software can fiuther comprise code to define two 
compartments, wherein one compartment includes biological processes related to a 
peripheral tissue environment and the second compartmmt includes biological processes 
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related to a lymphoid tissue environment. Further, the computer model can include a code to 
define the interaction between these two compartments. 

Additionally, the computer executable software code can comprise code to receive a 
user selection of a link representation &om a set of predefined link rq)resentations &om a set 

5 of predefined link representations, each predefined link representation being uniquely 
associated with a mathematical relationship. The user-selected link representation is 
associated with the interrelationship between the first biological variable and the second 
biological variable, a first link representation firom the set of predefined link represratations 
being a representation of ttie first biological variaLble having an effect on the second 

10 biological variable, a second link representation from the set of predefined link 

representations bemg a representation of instances of the first biological variable being 
converted to instances of the second biological variable. 

Another embodiment of the invention is a method for developing a compute model 
of the biological state of an adaptive immune response, comprising receiving user-selected 

IS . indications to define biological processes, each biological process being based on data that 
relates to changes in the adaptive immime response to biological attributes of a referaice 
pattern of adaptive immune response; producing a simulated biological attribute associated, 
with at least one biological attribute of the reference pattern of adaptive immune response; 
and assessing validity of tihie computer model based on a comparison between the simulated 

20 biological attribute and a corresponding biological attribute associated with the reference 
. pattern of adaptive immune response. 

Another, embodiment of the invention is a computer model of an ad^tive inunune 
response, comprising a computer-readable memory storing codes and a processor coupled to 
the computer-readable memory, the processor configured to execute the codes. The memory 

25 comprises code to define biological processes related to the biological state of the adaptive 
immune response and code to define math^atical relationships related to interactions 
among biological variables associated with the biological processes. At least two biological 
processes defined by the code are associated with the math^natical relationships. The 
combination of the codes stored in tihe memory that define the biological processes and the 

30 code that defines the mathematical relationships define a simulation of the biological state of 
flie adaptive immune response. 

The present invention also includes a method for developing an analytical model of 
an adi^tive immune response. This method includes the steps of identifying data related to 
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the biological state of the adaptive immune response; identifying biological processes related 
to the data, the biological processes defining at least one portion of the biological state of the 
adaptive immune response; and combining the biological processes to form an analytical 
representation of biological state of the adaptive immune response in the context of a 
peripheral tissue environment and a lymphoid tissue environment. 

Another embodiment of the invention is a method of developing an analytical model 
of the adaptive inmiune response wherein at least one biological process is associated with 
recruitment of immmie cells into the peripheral tissue environment. In a further embodiment, 
the method includes immune cells that are blood dendritic cells and blood monocytes. In yet 
a further embodiment, the method of developing a computer model of an adaptive immune 
response includes tiiie combining of biological processes so that the peripheral tissue 
environment is modeled with preferential recruitment of the blood dendritic cells over the 
blood monocytes. 

Jn one embodiment, in this analytical model, the analytical representation of the 
biological state of &e adq>tive immune response can be inq)lemented without the assistance 
of a computer systent 

Anoflier embodiment of the invention is a me&od for developing a computer model 
of an antigen*presCTting cell, conq)rising identifying data'relating to the physiological 

• regulatory mechanisms of the antigen-presenting cell, ttie data being associated wifli at least 
two fiom die group of antigen processing, migration, maturation, and mediator production of 
the antigCT-presenting cell and identifying biological processes related to the data, the 
biological processes defining at least one portion of tiie role of the antigen-presenting cell in 
an adaptive immune response. The biological processes are combined to form a simulation 

• offlieflmctioningoftheantigea-presentmgceUinthead£q)tiveimm bone 
embodiment, in the model the antigen-presentmg cell is a dendritic cell. In yet a furOier 
embodiment, the antigen-presenting cell is a myeloid dendritic cell.: In an additicmal 
embodiment, at least one of the biological processes is associated with a differential 
response to' antigen based on flie maturational state of the antigen-presenting cell. 

Components of the Mathematical Model 
Dendritic CeU (DC) Attributes 

DCs are one of several classes of professional antigen-presenting cells, which also 
include macrophages and B lynqphocytes. la one embodimrat of the invention, DCs are 
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modeled as the primary APC in a peripheral tissue. DCs continuously traf&c between the 
peripheral tissue and lymphoid tissue and function as sentinels of the immune system by 
finding and presenting antigens to lymphocytes. The DC mathematical model is designed 
to represent these biological attributes. 

In one embodiment, DC precursors are recruited into a peripheral tissue and 
differentiate into immature tissue DCs. The total population of DCs, comprised of DCs of 
all maturational states, is represented by a number of subpopulations which each have a 
discrete representative maturational state. The flux of DCs from the blood into a peripheral 
tissue is regulated by local environmental cues. Since the expression of adhesion molecules 
and chemoattractant receptors is dependent on maturational state, the rates of DC flux are 
roughly related to the rate of maturation. Once in the tissue, the flux of DCs between these 
diflFerent subpopulations and on to the LNs is regulated by the rate of maturation. 

The series of state nodes which represent the DC subpopulations in one embodiment 
is highlighted in FIG. 5 A and FIG. 5B. The interrelationship between nodes represents the 
Hux of DCs from the blood into the tissue and on to the LNs. In general, the maturation 
process, or the transition between subpopulation states, is characterized by distinct 
functional and cell surface matker changes, including downregulation of tissue homing 
chemokine recq)tors; antigen internalization, processing, and presentation capabilities; 
iq>regulation of costimulatory sur&ce molecules (e.g., CD80, CD86); and vpregulation of 
chemokine receptors for lynq)hoid tissue hommg. The late of DC maturation is modulated 
by the peripheral tissue microenviromnent, for example, tiuough cell-cell interactions (e.g. 
CD40 ligand-CD40 and Fas ligand-Fas) and through cytokines (e.g., BL-l, IL-3, II/4, IHO, 
TNF-alpha, and 6M-CSF). FIG. 6A and FIG. 6B higJiUght the areas of the DC modole in 
one embodiment that control the rate of DC influx to the peripheral tissue and maturation 
rate. 

The capability for an in^vidual DC subpopulation to ptocess antigen and activate T 
lymphocytes is dependent on the average maturational state of that subpopulation. 
lonnature DCs in the peripheral tissue are extremely efficient at internalizing and 
processing antigen, but are inefficient at stimulating T lyo^hocyte activatioa In contrast, 
mature DCs are efficient at initiating T lymphocyte activation, but are inefficient at 
intemalimg and processing antigen. The highlighted regions in FIG. 7 A and FIG. 7B 
represrat the combined effect of the antigen processing dependence on maturation and the • 
kinetics of antigen availability in the tissue. 
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Only DCs that have received the appiopriate stimuli are capable of activating 
antigen-specific T lymphocytes. Mature DCs express costimulatory molecules, including 
CD80 and CD86, and can potently activate T lymphocytes. Expression of costimulatory 
molecules is limited to mature DCs and is regulated by local environmental cues, such as 
the presence of IL-10. Mature DCs can present antigen to T lymphocytes in both peripheral 
tissues and lymphoid tissues. In a lymphoid tissue, DCs can present antigen, express 
costimulatory molecules, produce mediators that influence T helper lymphocyte 
differentiation, and be affected by T lymphocyte-DC cognate interactions. The DC ability to 
present antigen to T cells is dependent on the integrated effects of the kinetics of antigen 
availability in the tissue, ttie kinetics of DC maturation, and the ability of DCs to process 
antigen at discrete points during the maturational process. 

DC Cytokine Production 

DCs can produce mediators throughout tiieir lifecycle; this production is dependent 
on maturational state and local environmental cues. Mediators produced by DCs can 
include, for example, MDC, MCP-1, MlP-la, RANTES, TNF-a, DL^, E^S, II^IO, IL- 
12p40, and IH2p70. The nodes in FIG. 8 represent the regulation and production of DC 
mediators. The activity attributed to IL-12 is dependent on the concentration of the 
bioactive IL-12p70 heterodimer and the antagonists: II^12p40 homodimer and IL-12p40 
monomCT. The production of the possible IL-12 complexes (IL-I2p70 heterodimer, DL- 
12p40 homodimer, or IL-12p40 monomer) is dependent on the integrated exposure to local 
environmental cues during the DC maturational process. The highlighted regions of FIG. 
9 A , FIG. 9B, FIG. 9C and FIG. 9D represent this maturational dependence of the IL-12 
complexes. 

T Lymphocyte States 

' T lymphocytes and B lymphocytes are the two classes of cells diat possess 
* antigen-specific recognition receptors. The activation and subsequent expansion of 
antigen-specific cells subsequently lead to effector fimction in the immune response. In one 
embodiment of the invention, DC activation of CD4+ T he^er (Th) lymphocytes and their 
subsequent expansion and differentiation ^ere modeled. 

Presented with ttie appropriate antigen, in sufficient quantity, and qipropiiate 
costimulatory signals, antigeu-specific T lymphocytes can be activated to expand and 
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differentiate. The path of dijfferentiation may be regulated by the nature of the antigCT 
stimulus, cytokines^ and costimulatory molecules. Differentiated phenotypes have been 
identified and defined according to their pattran of cytokine production. By these 
definitions, a Thl population produces interferon-y (IFN-y) and interleukin-2 (E/-2), for 
example, while a Th2 population produces IL-4 and IL-5. Interestingly, data suggest that 
even under Thl or Th2 polarizing conditions, the T lymphocyte population may express 
multiple cytokines with different kinetics en route to becoming a polarized and largely 
^temiinally differentiated population. 

With this in mind, in one embodiment, the model can stratify the generic LN T 
lymphocyte population into subpopulations distinguished by characteristic patterns of 
cytokine productioiL Unpolarized, Thl, and Th2 populations exist, but these are subdivided 
to distinguish, for example, Th2 cells which produce IL-4 fiom Th2 cells which produce 
both IL-4 and IL-10. Different patterns of cytokine production by the T lymphocyte 
population as a whole can be realized by different distributions of cells among these 
subpopulations. 

FIG. 10 shows an exanq>le of an Bffect Diagram that displays a taxonomy of T 
lymphocyte states, according to an embodiment of the present invexxti^ Naive cells enter • 
and leave the LNs according to a steady state approximation. They can be made anergic or 
activated by antigen presentation. Following sufficient antigen stimulation, T lymphocytes 
produce IL-2 while in the primary effector I state (node labeled 'XN primary effector 1" in 
FIG. 10. Cells leaving the subsequent primary effector 2 state (node labeled primary 
fiffector:2'' in FIG. 10) will assume a Thl effector phenotype (node labeled *XN Thl 
effector^* in FIG. 10) or a Tli2 effector 1 phenotype (node labeled 'TN Th2 effector^ in FIG. 
10) according to the binding of their receptors for IFN-y and IL-4. Polarized Thl memory 
cells (node labeled Thl memory* in FIG. 10) and Th2 memory cells (node labeled 
•12^ Th2 memory** in FIG. 10) conqjete with naive cells for antigen presented on DCs, and 
are tiiemselves 8uscq>tible to anergy. A fraction of polarized effector T cells also traffic 
away fiom the LN, and a fraction of trafficking effector T cells subsequently gain entry to 
ithe peripheraltissue. 

Subpopulations of T cells are fiulher distinguished by varying suscq)tibilities to cell 
deatii in response to cytokine levels or cell-cell contact For example, the modeled Hil 
effector cells undergo activation-induced cell death (AICD) more rapidly than the modeled 
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Th2 effector cells. The model can include AICD and/or growth factor-withdrawal-induced 
apoptosis (GFWA). 



T Lymphocyte Cytokine Production 
5 Several cytokmes produced by LN T lymphocytes are involved in regulation of T 

, . lymphocyte, B lymphocyte, and DC LN populations. These can include, for example, IL-2, 
. IL-4, IL-5, IL-6, IL-10, and IFN-y. In addition, T lymphocytes that trafBc to a peripheral 
tissue produce cytokines and chemokines involved in the regulation of T lymphocytes and 
other cell populations (e.g., macrophages, eosinophils) in a peripheral tissue. These can 
10 " include, for example, IL-4, IL-5, 11^16, IFN-y, MIP-ip, and 1-309. FIG. 1 1 shows an 

example of a portion of an Effect Diagram that calculates T cell cytokine production. Badh 
defined T lymphocyte state produces a particular profile of cytokines at a specified rate. The 
kinetics of production are modified by the rate at Mf/iAch cells move fi:om one state to the 
next 

15 

T Lymphocyte Surface Molecule Ei^resslon and Ligation 

In obe embodiment, the modeled T lymphocytes can express several cell surface 
molecules involved in the regulation of various cell populations. The computer model can 
accommodate cell surface molecules involved in T lymphocyte regulatory Amotion 
20 including costimulation (e.g., CD28), cellular activation (e.g., CD40 ligand), and cell death 
(e.g., Fas ligand). 

T Lymphocyte Regulation of Expansion, Differentiation and Apoptosis 

. FIG. 12 shows an example of a portion of an Effect Diagram relating cytokine 
25 binding and cell-cell interactions to T lymphocyte expansion, differentiation, and death, 
according to an embodiment of the present invention. T lymphocyte expansion and 
differentiation are modified by cytokine growth factors, cytokine inhibitory factors, and 
• costimulatioiL T cell death is modified by insufficient concentration of cytokine growrth 
factors (e.g., IL-2, IL-4) or by activation-induced cell death (AICD). AICD occurs 
30 following ligation of death receptors (e.g., Fas) on the activated T cell surface. 

FIG. 13 shows Ihe Effect Diagram of FIG. 12 wifli an overlay that shows the 
equation used to calculate LN pe2 expansion function, according to an embodiment of the 
present invention. As shown in FIG. 13, an equation for this Amotion incorporates cytokine 
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effects, costimulatory molecule effects, cell cycle time, and a rate of movement for cells 
fix)m the pe2 state to ttie next state. 



T cell Trafficking 

Some fraction of the effector T lymphocyte populations traffic away from the LN. 
Some of the trafficking T lymphocytes will enter peripheral tissues where they are regulated 
and may affect other cell populations. T lymphocyte entry to the peripheral tissues is 
regulated by endoflielial expression of several adhesion molecules and by concentration 
gradients of chemotactic mediators. Cytokines induce endotheUal cell expression of several 
adhesion molecules. These may include P-selectin, E*selectin, vascular cell adhesion 
molecule-1 (VCAM-1), and intercellular adhesion molecule-1 (ICAM-1). Chemotactic 
mediators are produced by several cell types in the peripheral tissues and may include, for 
exampb, IL-8, IL-16, MIP-lp, thymus and activation regulated chemokine (TARC), and 
RANTES. Chemotactic mediators may have differential effects on Thl and Th2 entry to 
the peripheral tissues as Thl and Th2 cells express different sets of recq)tors. FIG. 14 
shows an Effect Diagram that depicts Thl and Th2 cell trafficking as modulated by 
di^otactic mediators and adhesion molecules. 

Emergent Regulation of Thl/Thl Polarization 

It has been shown that Thl/Th2 polarization can be influaiced by antigen dose, such 
that moderate levels of antigen produce a Thl response while high levels of antigen produce 
a 1112 response (Hoskenetal., J. Exp. Med. 182: 1579-1584, 1995). Similar data Mist 
.lioking sudi trends to increased costunulation and int^cellular adhesion. 

In one embodiment of ttie present invention, the computer model does not 
incorporate explicit signaling for naive cells to produce Thl or Th2 progeny with varying 
antigen levels. Nevertheless^ tins trend is observable as a result of differences in the rates at 
which Thl and Th2 cells e3q)and and undergo apoptosis followmg activation. For example, 
FIG. 15 illustra^tes an ou^ut of the model showing Thl and Th2 cells after primary culture 
as a fimction of antigen dose, according to an embodiment of the present invqition. As 
FJG. 15 shows, the number of Thl cells decreases with the antigen dose unlike the number 
of Th2 cells which increases with fixe antigen dose. 

In another embodiment, the levels of costimulation and - in a rudimentary way - 
intercellular adhesion can be varied within the conq)uter model to produce similar results. 
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Thus, explicit signaling is not required to observe these shifts in Thl/2 polarization with 
variations in these factors, and may be reinforced in nature by T cell population dynamics. 
Published data suggest that varying levels of antigen can produce multiple outcomes in T 
lymphocyte polarization under different experimental conditions. 

FIG. 16 illustrates an output of the model that depicts evolving numbers of Tbl and 
Th2 cells during a simulated primary in vitro culture at relatively low antigen dose, 
according to an embodiment of the present invention. Although the response shown 
ultimately &vors a Thl phenotype, Th2 cells predominate throughout early population 
development. As antigen dose is decreased, the phenotypic reversal of such a response 
occurs later and later in time as shown in FIG. 17. Depending on the point at which cell 
numbers are measured, an apparent Th2 response will therefore be obs^ed for sufficiently 
low levels of antigen. 

E3:amDle Development of a Model Component; Dendritic Cell Precnrsor Recruitment 

to a Peripheral Tissne 

The following discussion provides an txsmple of a method by which &e 
conq>oneDts of the above-described mathematical model can be developed. As discussed 
above, the various elements of the physiologic system are represented by the components 
shown in ihe Bffect Diagram. These conq>onmts are denoted by state and function nodes, 
which, with flie arrows and modifiers, represent mafiiematical relationships that define the 
elements of flie physiologic system. In general, these math^natical relationships are 
developed with the aid of pubhcly available or privately generated information on ttie 
relevant physiological components. The development of the mathemadcal relationships 
underlying the module diagram for DC precursor recruitment into the peripheral tissue will 
be discussed here as an example. 

FIG. 18 shows an example of an Effect Diagram for the recruitment of DC 
precursors into a peripheral tissue, where that peripheral tissue is the lun^. For illustration 
purposes, this Effect Diagram is a simplified version of the Effect Diagram for DC function 
depicted on page A-2 in Appendix A. The primary focus of this simplified mathematical 
model is to calculate the contributions of myeloid lung DC precursor populations, which 
exist in the blood, to the homeostatic and inflammatory response of the tissue DC 
populations in the human lung. The more detailed mathematical model depicted on page 
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page A-2 in Appendix A also includes the effects of additional mediators on the regulation 
of DC populations and additional states representing lung DCs. 

As FIG. 18 illustrates, maintenance and enhancement of lung DC populations in 
response to inflammatory stimuli are comprised of the regulated transport of three cell 
populations: node 1800, blood DC (BDC); node 1802, blood monocytes (BMo); and node 
1804, lung DC (LDC). The following discussion relates to deriving the underlying 
mathematical relationships for these physiological components based on the q)propriate 
publicly available information. Although not discussed herein, the remaining regulation of 
DC recruitment and population dynamics can be similarly derived from pubhcly available 
information. 

In the mammalian lung, DCs form an overlying network at the lung surface that has 
been associated with their role in acquiring and processing antigens at the air-lung surface 
interface. The total numb^ of LDCs changes dynamically when the lung is exposed to 
antigen; lung exposure to viial particles, bacteria, or allergen has hem shown to ten^rarily 
induce vp to 3.S-fold inoreases in the LDC density (McWiUiam, A. S. et al., JEjq? Med, 
179:1331-1336, 1994; Jahnsen, F. L, et al., ThoraXy 56:823-826, 2001). To represent fliese 
dynamics, the precursor populations for LDCs, the processes gbveming cellular influx into 
lung tissue, and the LDC numbers have been modeled. 

DC recruitment and associated population dynamics can be quantitatively 
represented as a set of coupled ordinary dififerential equations using basic engmeeriag 
principles such as a conservation of species (Bird, R. B. et al., Transport Phenomena^ 2"*^ 
Edition (2002), J. Wiley). In general, the relationships for the population dynamics of the 
three cellular classifications are: 

Rate of change of BDC = rate of BDC synthesis - rate of migration of BDC into 
lungs and other tissue compartments 

Rate of change of BMo = rate of BMo synthesis - rate of migration of BMo into 
hmgs and other tissue compartments 

Rate of change of LDC = rate of recruitment of BMo + rate of recruitment of BDC - 
rate of migration of LDC to lyn:q>h nodes. 
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lltese three ceU populations can be normalized to their total population sizes and 
then the rate of change of the populations can be written as the following equations: 

dBDC BDC 



dt BDCjoT 
dBMo _ BMogyf, 
dt BMoj„ 



-kfAM^it)'d'BDC-k,AMo-BDC (1) 



dLDC . 



- k, • AM I (/)•/?• BMo -k,AM'o -BMo (2) 



■m-.s-BDC+^^'fi'SM^-k^'WC (3) 

TOT ^^^TOT . ) 



.where 



AMj^{t) = time - dependent endothelial ligand in lung which facilitates recruitment (sites) 
AMq = endothelial ligand in other tissues which facilitates BDC recniitment (sites) 
AM^Q = endothelial ligand in oflier tissues which £icilitates BMo recruitment (sites) 
fi s= sdnsitivity of blood monocytes to endothelialligand expiessicm 
S = sensitivity of blood dendritic ceUs to endothelial ligand expression 
BDC = normalized blood dendritic cells (cells/total cells) 
BMo = normalized blood monocytes (cells/total cells) 



LDC = normalized lung dendritic cells (cells/total cells) 
.BDCsYff = synfliesis rate of blood dendritic cells (cells/hour) 
BDCj^ = total nmnb^ of blood dendritic cells (total cells) 
BMosm = synthesis rate of blood monocytes (cells/hour) 
fiMoj^Y = total number of blood monocytes (total cells) 
LDCj^ = total number of lung dendritic cells (total cells) 
kj = migration influx rate constant (l/(sites hour)) 
' kj^ = maturation rate constant (1/hom) 

Equations (1) and (2) define the dynamic state of BDC and BMo populations as 
regulated by rates of synthesis and loss to tissue compartments. Equation (3) defines the 
dynamic state of LDC^s as regulated by influx to and efQux fix>m the lung tissue. The first 
two terms on the rig}it-hand side of the equation (3) imply that once BMo aire recited into 
the tissue compartment they become LDCs. In one embodiment of ttie invention, the 
alternative pathway of monocyte differentiation into otfa^ cell types (i.e., macrophages) in 
the lung has not been represented The maturation or efflux rate of LDCs firom the lung 
compartment is the third term on tiie rigiht-hand side of equation (3). This rate of maturation 
of LDCs represents the migration of mature LDCs into the lymphatic system and. 
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ultimately, to the lymph nodes. This term also includes other pathways for the eUmination 
of LDCs such as apoptosis. 

The numerical output of these calculations determines the relative density of LDC, 
BDC, and BMo at any given time. The values of these variables can then be related to the 
uptake of antigen, antigen presentation, and other DC activities within the model. 

The limg endothelial adhesion molecule expression is dependent on the 
inflammatory state of the lung, hi the short term, the endothelial adhesion molecule 
expression in other tissues will not be affected by the lung inflammatory response. Thus the 
kinetics of recruitment to other tissues can be neglected. 

The parameters AM'o, kj, and ku can be eliminated using a steady-state 
{q>proximation for the BDC and BMo populations and equivalent cellular tumover rates. 
Equations (1) - (3) can then be expressed as: 

dBDC ^ BDCsm 
dt ADCto- 
dBMo BMosTf, 



dt BMOjoT 



dLDC _ A M 
dt 5 AM, 



fi P-AM,{t)H^-P)AM,it = 0) + AMo —\ 

^ [E^.S.BDC-.^^-fi-BM^ 

{t = 0) + AMoyLDC^- LDCjor ) 



AM^(t = 0) iBDCsrs S + BMo^.p) -—^ 
LDCj„i5AM^(t = 0) + AMo) 



(6) 



In general, parameter values can be defined, based on appropriate measurements in 
humans. However, in the absence of specific kinetic measurements in humans, data fix)m 
rodent or other studies were used. The estimates for parameter values for which no 
experimental data are available can be detemMued using by the solution of these ordinary 
differential equations in conjunction with experimental data that represents the appropriate 
general behaviors to be reproduced. Summed squared error of measurements can be used to 
determine goodness-of-fit 



ExpCTimentallv-defined parameters 

B:q>erimentally-defined paramet^ are ttiose whose values could be identified or 
derived directly fipom experimental data. These parameters are summarized in Table 1 and 
discussed below. 



34 



wo 03/001891 PCT/US02/20672 



Table 1. Experimentallv-defined parameter values 



Parameter Name 


Value 


BMotot 


2.8xl0^ceUs 


BDCtot 


9.744 X 10' ceUs 



The total population of blood monocytes (BMotot) and blood DC (BDCtot) were 
estimated to be 2.8x10^ cells and 9.744x10^ cells, respectively, based on published 
estimates of BMo density, BDC density, and assuming tiiat &e average human has 5.6 liters 
of blood (70 kg person, 8% of body weight is blood (Guyton, A.C., Textbook of Medical 
Physiology, 7* Ed., 1986, W.B. Saunders Company), and density of blood is 1 g/ml). The 
published estimate of BMo density is S.QxlO^ cells/L (Rich, L N., Monocytes and 
Macrophages in Primary Hematopoietic Cells (1999) Kluwer Academic Publishers) while 
that of BDC density is 17-4±5.4xl0* cells/L (Upham, J. W. et al.. Cytometry, 40:50-59, 
2000). 

Driving function which represents mf lflnwifltnry gtinmli 

By executing the mafhematical model with a driving function, AML(t), which 
reflects a dynamic inflammatory stimulus, flie equations may be used to determine dynamic 
changes in recruitment of precursor D C populations into the lung tissue and to determine 
dynamic increases in the LDC population under inflammatory conditions. The specification 
of AML(t) is explained here. Based on the rapid recruitment of tiie blood cell populations 
into the lung, P-selectin is the adhesiori molecule with expression kinetics most capable of 
iiiitiating recruitment. P-selectin is stored in Weibel-Palade bodies in endothelial cells and 
has been observed to increase several-fold quickly on the cell surface, reaching a maximum 
. after 5-10 minutes (e.g., Tedesco, F. et al., J Exp Med, 185:1619-1627, 1997). 
Disappearance of the surfiice P-selectin via «idocytosis has been observed within 30-60 
minutes. In addition, P- and B-selectin can both be transcriptionally induced by 
inflammatory mediators with similar kinetics and peak approximately 14 hours following 
stimulation (Vestweb^, D. and Blanks, J. E., Physiol Rev, 79:181-213, 1999). To capture 
these kinetic features, the driving fimction AML(t) has been estimated to have the time- 
course shown in FIG. 19. The value of AML(t) at time equal to zero (AML(t=0)) was set to 
0.05683 sites (Table 2). 
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Table 2. Estimated parameter values for mathematical model 



Parameter Name 


Value 


AML(t = 0) 


0.05683 sites 



Simulation-defmed parameters 

Simulation-defined parameters are those whose values were selected to yield 
simulation results that were substantially consistent with experimentally known behaviors 
of DC precursors in the blood and LDCs. The process of selecting these parameter values 
involves reproducing experimratal protocols in the mattiematical model and optimizmg 
parameter values based on the goodness-of-fit of simulation results to expramental results. 
The simulation-defined parameters are listed in Table 3 and described below. 

Table 3* Simulation-defined parameter values. 



Parameter Name 


Value 


BMOSYN 


2.526 X 10' cells/hour 


BDCsYN 


8.792 X 10' cells/hour 


AMo 


0.02632 sites 




1 


5 


. 43.07 


XDCtot 


2.739x10' cells 



Using data jGx>m published ladiolabeling Studies (Wluteto^ 
464, 1966), the turnover rate of BMos was estimated to be 76.8 hours. The experimental 
data for this radiolabeling study are shown in FIG. 20, and corresponding simulation results 
using tills half-life are also shown. Since both the BMo and BDC come fix)m a common 
precursor population, the BDC turnover rate was assumed to also be 76.8 hours. 

Using a steady-state approximation with equation (1), and this turnover rate, the 
synthesis rate of blood monocytes (BMosyn) was calculated to be 2.526x10^ cells/hour. 
The same process with equation (2) gives a calculated synthesis rate of blood DCs 
(BDCsyn) of 8.792x10^ cells/hour. 

Values for AMq, S, and LDCtot were determined fix>m data on steady-state and 
inflammatory chall^ge studies of LDCs (Upham, J. W. et al.. Am JRespir Crit Care Med, 
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159:A854, 1999; McWiUiam, A. S. et al., J Exp Med, 1994, 179:1331-1336; Holt, P. G. et 
al., JImmtmoly 1994, 153:256-261). The accuracy of the simulation-defined parameters is 
evaluated against the ability of the mathematical model to reproduce the experimental data. 
The data used to set the parameter values include the following: Upham et al. reported that 
5 BDC and BMo dropped approximately 37% and 5%, respectively, at three hours following 
antigen challenge. After 24 hoiurs, the BDC level had returned to 78% of the pre-challenge 
level while the BMo level had returned to normal. McWilliam et al. observed IDC numbers 
peak at 3.5x.baseline values within 1 hour after introduction of an inflammatory stimulus. 
' Holt et al. demonstrated the decline of LDC populations following irradiation of bone 
10 maCTOw. 

. . FIG. 2 1 demonstrates that using the BDC and BMo values defined in Tables 1 -3 and 
equations 4-6 with a simulation protocol that mimics the experimental protocol of Holt et 
al., the LDC values are in agreement with the e}qperimental findings of Holt et al. In FIG. 
22, title simulation results for both the blood and tissue populations are shown to be 
IS substan^ally consistent with the experimental data of Upham et al. and McWilliam et al. 
The results are also substantially consistent with other measurements of the dynamic 
response of the LDC populations to various experimental conditions (Ibgulli, E. et aL, J&qf 
Med, 1997, 185:2133-2141; Lambiecht, B. N. et al., Am JRespir Cell Mol Biol, 1999, 
20:1165-1174; Stumbles, P. A. et vl., J Immunol 2001, 167:228-234). 

20 

Useful outcomes of mafliCTiatical modeliiig 

Malhematical modeling of DC recruitment as just explained has generated insightful 
observations on this area of biology. Upham et al. demonstrated that at 3 hours following 
antigen challCTge, BDC density dropped by 37% while BMo density dropped by 5%. In 

25 rq>roducing these experimmtal data with the model, it was also found that BDC account for 
20% of the post-challenge tissue DC population. So, while BDC comprise only about 3.4% 
of the total DC precursors in the blood {i.e., BDC plus BMo), they comprise a much larger 
percentage of the lung tissue population. The functional consequences of preferential BDC 
reaiiitment may be profound. BDC and BMo are functionally distmct in theu: antigen 

30 processing and presentation abilities (Caux, C. et al.. Blood, 90: 1458-1470, 1997; Gairett, 
W. S. et al.. Cell, 102:325-334, 2000; Sallusto, R and Lanzavecchia, A., JE)q> Med, 
179:1 109-1 1 18, 1994; Yang, D. et al., J Immunol, 165:2694-2702, 2000). Hence, 
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preferential recruitment of one cell type over the other may dramatically influence antigen 

presentation to T lymphocytes. 

An additional useful observation is that the kinetics of DC precursor and tissue DC 

population responses following antigen challenge suggest that P-selectin is the most likely 
5 . among the known adhesion molecules to mediate recruitment because it can be rapidly 

mobilized to the endotheUal cell surface upon stimulation (Tedesco, F. et al., J Exp Med, 

185:1619-1627, 1997). 

As this example model of DC precursor recruitment into a peripheral tissue 

generally illustrates, the components of the Effect Diagram, denoted by state and function 
10 nodes and the arrows and modifiers linking them, represent mathematical relationships that 

define the elements of the physiologic system. These mafliematical relationships can be 

developed with the aid of appropriate publicly available information on the relevant 

physiological components, hi other words, the Effect Diagram indicates the types of 

mathematical relationships tibat are represmted in tiie model. The publicly available 
15 . information can be put into a form that ntiatc^es the stnicture of die Bffw^ Ihdiis 

way, the stiiicture of fhe model can be developed. 

Simulation of Biological Attributes of an Adaptive Immune Response 
The following discussion describes the nature of the biological attributes that can be 
20 obtained by numerical or analytical integration of the mathematical model. It further 

elucidates changes that may be made to fhe model to obtain simulated biological attributes 
that correspond to quahtatively or quantitatively different ad^tive inunime responses. 

The mathematical model is equipped with a set of baseline parameters selected to 
represoit a particular type of adaptive immune response. In one embodiment, tfie baseline 
25 parameters are selected such that the simulated biological attributes are substantially 

consistent with the biological attributes of an established immune response to an allergic 
stimulus. The parameters of the model can be changed to represent varying manifestations 
of fhe response including for example, acute responses to a bolus e^qposure of antigen, or 
low level chronic responses to low levels of antigen, or quiescent response to fhe absence of 
30 antigen. 

The model can also be changed parametrically to represent differrat contributions of 
fhe involved biological processes to fhe biological state. Changing fhe contributiQns to the 
biological processes will yield different simulated biological attributes and enable 
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exploration of how parameter changes affect outcomes. For example, changing the 
appropriate model parameters such that DC bioactive IL-12 production is enhanced favors 
Thl polarizing conditions and enables investigation of the subsequent changes imposed on 
T lymphocyte populations. 

When integrated with other components to generate a normal or diseased 
physiological system (e.g., allergic asthma), the model can represent the contribution of the 
immune response to the state of the system or to the progression between disease states of 
different severity. The model can also represent the different stages of immune responses, 
i.e., primary and secondary exposure to a given antigen, as well as exposure to varying 
levels of antigen. One means of generating this variation in the model can involve rq)lacing 
one or more biological variables, fomoierly fixed at a particular value, with one or moie 
•biological variables fliat evolve over time and depend on some previously included or new 
biological processes. For example, in one embodiment, the APCs are exposed to a chronic 
low level of antigenic stimuli. The low level of antigen exposure results in some basal level 
of inflammation in flie peripheral tissue and maintains a chronic low level immune * 
response. Altering the immune response, to include for example, acute responses, might 
involve Tq>lacing a fixed parameter (e.g., antigen level) with a direct fiinction of time, an 
algebraic function of other biological variables (i.e., a biological process), or via a 
dynamical systems equation such an ordinary differential equation. Alternatively, changing 
tfie immune response might involve adding new processes such as exposure to viral or 
bacterial pathogens, influx and dynamics of NK cells and neutrophils, and allowing APC 
behavior to be modulated by these processes. 

In one embodiment the previously fixed values that specify allergen exposure of a 
peripheral lung tissue by a mild allergic asthmatic, are rq>laced by a direct function of time 
or by a function of other biological variables to represent the effect of seasonal changes in 
allergic stimuli to the pathology of allergic asthma. The model depiction of acute 
exacerbations in a chronically inflamed lung can be used to study, for example, the role of 
the adaptive immune response in acute exacerbations as well as chronic inflammation and 
q)proaches to alter the characta of the ad^^tive inmoiune response through ther^eutic 
changes to the APCs themselves or the peripheral tissue environment 
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Model Calibration 

Model calibration refers to the estimation of parameter values based on quantitative 
and qualitative experimental observations tiiat cotrespond to biological variables, biological 
processes and/or biological states represented in the mathCTiatical model. Parameter values 
are estimated from both in vitro and in vivo data from the Uterature in the model calibration 
process. With parameter values estimated in this way, the modules of the mathematical 
model show behavior substantially consistent with experimental studies that focus on 
specific aspects of the adq)tive immune response. For example, the model can reproduce 
the timing and numbers of cells involved in DC recruitment into peripheral tissue, and the 
response of different T lymphocyte populations to APCs with varying antigen loads. As 
shown in the previous section, '"Example Development of a Model Component: Dendritic 
Cell Precursor Recruitment to a Peripheral Tissue", the recruitment of DCs has been shown 
to match both blood and periph^ tissue data. Published data specifying the q)piopriate 
behavior of the system overall can be used to calibrate remaining degrees of freedom in the 
computer model. 

Examples of tiie calibration of the migration kinetics of DCs into the lymph node are 
shown in FIG. 23 and FIG. 24. As can be seen in the FIG. 23, the simulation result is 
substantially consistent with experimental measurements reported by two independent 
groups in the published literature (lagulli et al., JEjq? Med, 185: 2133-2141, 1997; 
Vermaelen et al., J Exp Med, 193: 51-60, 2001) in the context of non-productive 
interactions with CD4+ T lymphocytes, hi contrast, productive interactions of DCs with 
CD4+ T cells result in sUghtly different kinetics, as seoa in FIG 24. Productive interactions 
induce activation of the CD4+ T lymphocytes, which in turn increases the q)optosis of the 
antigen-presenting DCs. The simulation result is consistent with pubUshed experimental 
results (Ingulli et al., J Exp Med, 185: 2133-2141, 1997) and is shown in FIG. 24. 

The stimulation of T lymphocytes by APCs also reproduces the experimental data 
reported in the scientific litCTature. In Table 4, the simulation resxilts are con:q)ared to a 
study by London et al. (London et al., JImmmiol, 164:265-272, 2000) wha-e the cytokine 
production response of various T lymphocyte populations to antigen and costimulation was 
measured. The model parameters were selected to be consistent with the experimental 
protocol described by London et al. in which both memory (mem) and naive T lymphocyte 
populations were stimulated with two levels of antigen (Ag) and costimulatory molecules 
expressed by APCs. The cytokine production kinetics were measured (Table 4), and in 
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some instances several cytokines were not observed above the detection limit (n.d.). The 
simulation results demonstrate that both naive and memory populations in the mathematical 
model respond appropriately in terms of the relative cytokine amount and timing to different 
levels of antigen. The general behavioral diiSerences between naive and memory T 
lymphocyte responses are also consistent with the understanding that memory lymphocytes 
require less antigen for stimulation as demonstrated by the data of London et al. 



TABLE 4. 



Cytokine 


T Cell phenotype and 
stimulation 


Simulation Results 


Results from London et 
al. 2000 






Nonnalize 

dpeak 

value 


Time of 
peak(h) 


Nonnalized 
peak value 


Time of 
peak(h) 


TL-2 


Na^re, low Ag 


<0.001 




n.d. 




JL-2 


Narvre,liighAg 


0.053 


62 


0.04 


60 


TL-2 


Mem, low Ag 


0.122 


58 


0.07 


60 


JL-2 


Mem,Iiig^Ag 


0.157 


56 


0.16 


60 


IL-4 


NalVe (Mgh and low) 


<0.001 




n.d. 




IL-4 


Mem, low Ag 


0.007 


60 


0.005 


60 


IL-4 


Mem, high Ag 


0.008 


58 


0.008 


84 


IFN-ganuna 


Naive (high and low) 


<0.001 




ad. 




IFN-gamma 


Mem, low Ag 


0.015 


64 


0.018 


60 


n^-gamma 


Mem, high Ag 


0.018 


70 


0.018 


84 



Tnitialfaation^f the Mathematical Eaaations and Numerical Solution of the Computer 

Model 

Since the Effect Diagram defines a set of ordinary diffecential equations as described 
above, once the initial values of the biological variables are specified, along with the values 
for the model paramet^, the equations can be solved niunerically by a computer using 
standard algorithms. See, for example, William H. Press et al. Numerical Recipes in C: The 
Alt of Scientific Computing, 2nd edition (January 1993) Cambridge Univ. Press. As 
illustrated above in the Example Development of a Model Component: Dendritic Cell 
Precursor Recruitment to a Peripheral Tissue section, one can derive equations, obtain 
initial conditions, and estimate parameter values fix)m the public literature. likewise, other 
initial conditions and parameter values can be estimated for different conditions and can be 
used to simulate the chronological progression of the biological state. 
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In one embodiment, the computer executable software code numerically solves the 
mathematical equations of the model under simulated experimental conditions. For 
example, one could simulate an in vitro experiment by specifying the duration of the 
experiment and the following initial conditions for the biological variables: DC density, the 
DC state(s), the T lymphocyte density, the T lymphocyte population states (e.g., naive 
GD4+ T lymphocytes, resting Thl lymphocytes, resting Th2 lymphocytes), the amoimt of 
antigen, and the cytokine environment (e.g., endogenous or exogenous IL-12, IFN-gamma, 
JLr4), The numerical solution would include tiie values for all the experimentally measured 
cell populations and mediator levels (e.g., number of ILr2 expressing Thl lynq)hocytes) at 
the times' they were measured in the laboratory. In addition, the nimierical solution could 
generate the complete chronological progression of all biological variables in the model 
over the course of the experiment. 

Furthermore, ttie computer executable software code can facilitate visualization and 
manipulation of the model equations and their associated parameters to simulate different 
patients subject to a variety of stimuli. See, e.&, U.S- Patent 6,078,739, entitled "Managing 
objects and parameter values associated wi& ttie objects wi^ a simulation modeV' and 
U.S. Patent 6,069,629, entitled **Method of providing access to object parameters wifhin a , 
simulation model" the disclosures of which are incorporated herein by reference. 

In one embodiment the invention can be used to model therapeutic agents such as 
steroids, P-agonists, or leukotriene antagonists. Thus, the model can be used to rapidly test 
hypotheses and investigate potential drug targets or therapeutic strategies. For example, a 
therapy can be modeled in a static manner by modifying the Parameter Set of the 
q)propriate tissue(s) to represent the affect of the treatment on that tissue(s). Alternatively, 
tibierapeutic treatments can be modeled in a dynamic manner by allowing the user to specify 
the delivery of a treatment(s), for example, in a time-varying (and/or periodic) manner. To 
do this, the computer model has the ability to include pharmacokinetic representations of 
various modulatory classes of treatment (e.g., anti-cytokine antibodies, adjuvant-like 
mediators, steroids) and how these treatments can interact with the various cell types in a 
dynandc manner. Fuither, whea the model is integrated with a disease model, fliere is an 
ability to include pharmacokinetic representations of various tti«:sq>eutic classes (e.g., 
anti-cytokine antibodies, altered forms of antigen, adjuvant-like mediators, steroids) and 
evaluate how tixese fhen^eutics interact with tiie elements of peripheral and lymphoid tissue 
to generate a clinical outcome. 
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Validation of the Model 

The behavior of the model is validated by comparing simulated biological processes 
to reference patterns of those biological processes. Validation in this manner can be used to 
validate both the computer model and the analytical model embodiments of the invention. 
For example, one method of validating the behavior of the computer model is validated by 
comparing simulated biological attributes of the model to reference patterns of individual 
components of the adaptive immune response or to reference pattems of the entire adaptive 
immune response. An alternate method is to link the mathematical model of the adaptive 
immune response to a model of another biological syst^ that interacts with the adaptive 
immime response, for example a model of an allergic asthmatic limg, a model of intestinal 
bowel disease (IBD), or a model of Schistomsoma mansoni infection, and con^)are . 
simulations of the combined model to reference pattems for the combined system. 

In one embodiment, validation of the model of the ad^tive immune response is 
petfimned by linking it to a model of a peripheral tissue whose function depends on the 
ad^tive immune response, and compaiing simulation results for the linked model to 
teference pattems of the peripheral tissue biological function. The model of the adaptive 
immune response could also be integrated with any number of other conq)onents to 
represent aliealtfay or diseased physiological system. Validation of tiie tnodel would be 
done by comparing simulated results of the linked model to reference patterns of the 
components, which could be tissues or entire organisms that are dependent on the adaptive 
inunune response. 

To use linking of the adaptive immune response to a peripheral tissue model for 
. validation, the following process can be used The method of validation can involve having 
a model with the attributes of a peripheral tissue of interest in a particular biological state. 
Specific characteristics of cell types, cellular abnormalities, physiological abnormalities, 
and the chemical mediator environment of the peripheral tissue can be modeled 
qjpropriately. In particular, the peripheral tissue modeling is necessary to reflect the impact 
of peripheral tissue constituents on cells of the adaptive immune system and conversely, the 
inopact of immune cells on pmph^ tissue constituents. In one embodiment, the ad^tive 
immime response model can modified to reflect the nature of APCs and lymphocytes that 
are associated with specific biological state of flie selected peripheral tissue. The 
interactions of the ijamiune cells with each other and with cells or chemical mediators of the 
peripheral tissue result in simulated biological attributes that can be conq)ared to 
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experimentally observable reference patterns. Methods for validation of computer models 
are described in an application entitled "Apparatus and Methods for Validating a Computer 
Model," filed on May 16, 2002, Application Number 10/151,581 which is incorporated 
herein by reference. 

5 The adaptive immune response is imphcated in many diseases including allergic 

asthma, rheumatoid arthritis, inflammatory bowel disease, cancer, and all infectious 
diseases. As an example of validation by linking with another model, the adaptive immune 
response model can be validated througih incorporation in a disease model of allergic 
asthma. Under chronic conditions, the adaptive immune response model can, for example, 

10 provide the proper stimuli to produce known reference patterns in asthmatic patients 
including elevated IgB levels, partial degranulation of mast cells, airway 
hyperresponsiveness, and elevated levels of cytokines associated with a Th2 adaptive 
immune response. The model can also reproduce appropriate reference patterns of patient 
responses to bolus doses of antigen, including conqiiomised airway fimction and elevated 

15 . cell and chemical mediator levels in flie airways. 

. Consistency with module-specific reference patterns measured in in vivo and ' in vitro 
studies, as well as reference patterns of clinical outcomes when incorporated within a fiiU 
disease model, provides validation for die computer model. The mathematical model of an 
adaptive immune response can be considered a valid model if die simulated biological 

20 • attribute obtained is substantially consistent with a corresponding biological attribute 

obtained fix)m a cellular or whole animal model of an adaptive immune response. As the 
understanding of the adaptive immune response, and the diseases associated with the 
adaptive immune response, evolve in the art, the responses against which the model is 
validated can be modified. 

25 While various embodiments of the invention have been described above, it should be 

understood that they have been presented by way of example only, and not limitation; 
Thus, the breadth and scope of the present invention should not be limited by any of the 
• above-described ^bodiments, but should be defined only in accordance with the following 
claims and their equivalents. 

30 The previous description of the embodiments is provided to enable any person 

dolled in the art to make or use the invention. While the invention has been particularly 
shown and desqibed with reference to embodiments thereof it will be understood by those 
skilled in the art that various changes in form and details may be made tiierein without 
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departing from the spirit and scope of the invention. For example, although a certain 
embodiment of a computer system is described above, other embodiments are possible. 
Such computer system embodiments can be, for example, a networked or distributed 
computer system. 
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1 . A method for developing a computer model of an ad^tive immune response, 
comprising: 

identifying data relating to a biological state of the adaptive immune response; 

identifying a plurality of biological processes related to the data, the plurality of 
biological processes defining at least one portion of the biological state of the adaptive 
immune response; and 

combining the plurality of biological processes to form a simulation of ttie biological 
state of the adaptive immune response. 

2. The method of claim 1 , wherein the biological state of the adaptive immune 
response is a biological state of an acute response. 

3. The method of claim 1, wherein the biological state of the adaptive immune 
response is a biological state of a chronic response. 

4. The method of claim 1, wherein at least one biological process from the plurality of 
biological processes is associated with a biological variable that is a therapeutic agent 

5. The method of claim 1, further comprising: 

producing a simulated biological attribute associated wifli the biological state of the 
adaptive immune response; 

comparing the simulated biological attribute with a corresponding biological 
attribute in a reference pattern of the adaptive immune response; and 

identifying ttie con[tput^ model as a valid conoputer model of the adaptive immune 
response if the simulated biological attribute is substantially consistent with the biological 
attribute associated with the reference pattern of the adaptive immune response. 

6. The method of claim 1, wherein the combining the plurality of biological processes 
includes: 

forming a ISrst mathematical relation among biological variables associated with a 
first biological process from the plurality of biological processes; and 
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forming a second mathematical relation among biological variables associated with 
the first biological process and biological variables associated with a second biological 
process fix)m the plurality of biological processes. 

7. The method of claim 6, further comprising: 

5 creating a set of parametric changes in the first mathematical relation and the second 

mathematical relation; and 

producing a simulated biological attribute based on at least one parametric change 
from the set of parametric changes, the simulated biological attribute being substantially 
consistent with at least one biological attribute associated with a reference pattem of tibie. 

1 0 adaptive immune response. 

8. The method of claim 6, further comprising: 

converting a first biological variable into a converted biological variable a value of 
vMch changes over time, the first biological variable being associated with at least one 
fiom ike first matheinatical relation and the second mathematical relation; and 

IS . producing a series ofsimulated biological attributes based on the converted 

biological variable, the series of simulated biological attributes being substantially 
consistent with a corresponding biological attribute associated with a reference pattem of 
tiie adaptive immune response, the series ofsimulated biolo^cal attributes representing the 
chronological progression of the corresponding biological attribute in the reference pattem 

20 of the adaptive immune response. 

9. The method of claim 6, further conxprising: 

converting a parameter into a new biological variable a value of which changes over 
..time, the parameter being associated with at least one fix)m the first mathemadcal relation 
and the second mathematical relation; and 
25 producing a series of simulated biological attributes based on the new biological 

variable, the series ofsimulated biological attributes being substantially consistent with a 
biological attribute associated with a reference pattem of an ad2q)tive immune response, fiie 
series ofsimulated biological attributes representing the chronological progression of 
corresponding biological attributes m the reference pattem of the ad^tive immune 
30 response. 
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10. A mefbod for developing a computer model of an adaptive immune response, 
compiising: 

identifying data relating to a biological state of the adaptive immune response; 

identifying a plurality of biological processes related to the data, the plurality of 
5 biological processes defining at least one portion of the biological state of the adaptive 
immune response; and 

combining the plurality of biological processes to form a simulation of the biological 
state of the ad^^tive immune response in the context of a peripheral tissue environment and 
a lymphoid tissue environment. 

10 11. The method of claim 10, wherein at least one biological process fiom the plurality of 
biological processes is associated with recruitment of immune cells into the peripheral 
tissue environment. 

12. The method of claim 1 1 , wherein the immune cells are blood dendritic cells and 
blood monocytes. 

15 13. The method of claim 12, wherein the plurality of biological processes are combined 
so that the perqpheral tissue aivironment is modeled with preferaitial recraitment of the 
blood dmdritic cells over the blood monocytes. 

14. The method of claim 10, wherein at least one biological process j&om the plurality of 
biological processes is associated wifli a biological variable that is a tfier^eutic agent. 

20 

15. The me&odofclaim 10, further conqncising: 

producing a simulated biological attribute associated with the biological state of the 
adaptive immune response; 

comparing the simulated biological attribute with a corresponding biological 
25 attribute in a reference pattern of tiie adq>tive immune response; and 

idmtifying ttie computer model as a valid computer model of the adaptive immune 
response if the simulated biological attribute is substantially consistent with the biological 
attribute associated with the reference pattern of the adaptive immune response. 



16. The mettiod of claim 10, wherein the combining the plurality of biological processes 
30 includes: 
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fonning a first mathematical relation among biological variables associated with a 
first biological process from the plurality of biological processes; and 

fonning a second mathematical relation among biological variables associated with 
the first biological process and biological variables associated with a second biological 
process from the plurality of biological processes. 

17. The method of claim 16, further comprising: 

creating a set of parametric changes in the first mathematical relation and the second 
mathematical relation; and 

producing a simulated biological attribute based on at least one parametric change 
fix)m the set of parametric changes, the simulated biological attribute being substantially = 
consistent with at least one biological attribute associated with a refermce pattern of the 
adaptive immune response. 

18. The method of claim 16, further comprising: 

converting a first biological variable into a converted biological variable a value of 
which changes over time, the first biological variable being associated with at least one 
from the first mathematical relation and the second mathematical relation; and 

producing a series of simulated biological attributes based on the converted 
biological variable, the series of simulated biological attributes being substantially 
consistent with a corresponding biological attribute associated with a reference pattem of 
the adaptive immune response, the series of simulated biological attributes representing the 
chronological progression of the corresponding biological attribute in the reference pattem 
of the adaptive immune response. 

1 9. The method of claun 16, fimher comprising: 

converting a parameter into a new biological variable a value of which changes over 
time, the parameter being associated with at least one firom the first mathematical relation 
and the second mathematical relation; and 

producing a series of simulated biological attributes based on the new biological 
variable, the series of simulated biological attributes being substantially consistent with a 
biological attribute associated with a reference pattern of an adaptive immune response, tfie 
series of simulated biological attributes representing the chronological progression of 
corresponding biological attributes in the reference pattem of the adaptive immune 
response. 
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20. A computer model of an adaptive immune response, comprising: 

code to define a set of biological processes related to a biological state of the 
adsqitive immune response; and 

code to define a set of mathematical relationships related to interactions among 
biological variables associated with the set of biological processes, at least two biological 
processes from the set of biological processes being associated with the set of mathematical 
relationships, a combination of the code to define the set of biological processes and the 
code to define the set of mathematical relationships defining a simulation of the biological 
state of the adaptive immune response in the context of a peripheral tissue environment and 
a lymphoid tissue enviroxmient 

21 . The computer model of claim 20, wherein at least one biological process from the 
set of biological processes is associated with reouitment of immune cells into flie peripheral 
tissue environment 

22. The computer model of claim 21, wherein ttie immune cells are blood dendritic cells 
and blood monocytes. 

23. The computer model of claim 22, wherein the set of biological processes are 
combined so that the peripheral tissue environment is modeled with preferential recruitment 
of the blood dendritic cells over the blood monocytes. 

24. The conq>uter model of claim 20, wherein at least one biological process fiom the 
set of biological processes is associated with a biological variable that is a therapeutic ag&it 

25. The computer model of claim 20, wherein upon execution of the code, a simulated 
biological attribute for the adaptive immune response is produced, the simulated biological 
attribute being substantially consistent with at least one biological attribute associated with 
a reference pattern of ttie adaptive immune response. 

lis. The computer model of claim 20, fiirther comprising: 

code to define a first compartment, said first compartment including biological 
processes related to aperipheral tissue environment, and 
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code to define a second compartment, said second compartment including biological 
processes related to a lymphoid tissue environment. 

27. The compute model of claim 26, further comprising: 

a code to define a set of biological processes related to trafficking of immune cells 
between said first and second compartments. 

28. A computer executable software code, comprising: 

code to define a plurality of biological processes related to a biological state of an 
adaptive immune response including: 

code to define a set of mathematical relations associated with a first biological 
process fix)m the plurality of biological processes and associated with interactions among 
biological variables associated vnih the first biological process, and 

code to define a set of mattiematical relations associated with a second biological 
process fix>m the plurality of biological processes and associated with interactions among 
biological variables associated wifii the second biological process, the plurality of biological 
processes being associated with the adaptive immune response in the context of a peripheral 
tissue environment and a lymphoid tissue environment. 

29. The computer executable software code of claim 28, wherein at least one biological 
process fi:bm tiie plurahty of biological processes is associated with recruitment of inmune 
cells into the peripheral tissue environmmt 

30. The computer executable software code of claim 29, wherein the immune cells are 
blood dendritic cells and blood monocytes. 

3 1 . The computer executable software code of claim 30, wherein the plurality of 
biological processes are combined so that the peripheral tissue environm^t is modeled with 
preierential recruitment of the blood dendritic cells over the blood monocytes. 

32. The computer executable sofhvare code of claim 28, wh^ein at least one biological 
process firom the plurality of biological processes is associated with a biological variable 
tiiat is a therapeutic agent 
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33 . The compute executable software code of claim 28, further comprising: 

code to receive a user selection of a link representation jfrom a set of predefined link 
representations, each predefined link representation in the set of predefined link 
representations being uniquely associated with a mathematical relationship from the set of 
mathematical relationships, the user-selected link representation being associated with the 
interrelationship between a first biological variable and a second biological variable, 

a first link representation fi'om the set of predefined link representations being a 
represen^tion pf the first biological variable having an effect on the second biological 
variable, 

a second link representation from the set of predefined link representations being a 
representation of instances of the first biological variable being converted to instances of the 
second biological variable. 

34. The computer executable software code of claim 28, further comprising: 
code to define a first compartment, said first compartment including biological 

processes relatied to a peripheral tissue environment, and 

' code to define a second compartment, said second compartmrat including biological 
processes related to a lymphoid tissue environment 

35. The computer executable software code of claim 34, further comprising: 

a code to define a set of biological processes related to trafficking of immune cells 
between said first and second compartm^ts. 

36. A method for developmg a computer model of an adaptive immune response, 
comprising: 

receiving a plurality of user-selected indications to define a plurality of biological 
processes, each biological process from the plurality of biological processes being based on 
data that relates changes in a biological state of the adaptive immune response to biological 
attributes of a reference pattern of the adaptive immune response; 

producing a simulated biological attribute associated with at least one biological 
attribute of the reference pattern of the adaptive immune response based on the combined 
plurality of biology processes; and 

assessing validity of the conq)uter model based on a comparison between the 

simulated biological attribute and a corresponding biological attribute associated with the 

reference pattern of the adaptive immune response. 

69 



wo 03/001891 PCT/US02/20672 

37. A computer system of an adaptive immune response, comprising: 
a computer-readable memory storing: 

code to define a set of biological processes related to a biological state of the 
adaptive immune response; and 

code to define a set of mathematical relationships related to interactions among 
biological variables associated with the biological processes, at least two biological 
processes from the set of biological processes being associated with the set of math^atical 
relationships, a combination of the code to define the set of biological processes and the 
code to define the set of mathematical relationships defining a simulation of the ads^tive 
immune response in the context of a peripheral tissue environment and a lymphoid tissue 
environment; and 

a processor coupled to the computer-readable memory, the processor being 
configured to execute the codes. 

38. The computer system of claim 37, wherein at least one biological process horn the 
set of biological processes is associated with a biological variable that is a therapeutic agent. 

39. The computer system of claim 37, wherein upon ^ecution of the code, a simulated 
biological attribute for flie ad^tive immune response is produced, the simulated biological 
attribute being substantially consistent with at least one biological attribute associated with 
a reference pattern of the adq)tive immxme response. 

40. The computer system of clahn 37, fiirther comprising: 

code to define a first compartment, said first con9)artment includes biological 
processes related to a p^pheral tissue environment; and 

code to define a second compartment, said second compartment includes biological 
processes related to a lymphoid tissue envirpnment 

41. The computer system of claim 38, further comprising: 

a code to define a set of biological processes related to trafficking of immune cells 
between said &st and second compartments. 

42. A metfxod for developing an analytical model of an adqitive immune response, 
comprising: 

identifying data relating to a biological state of the ad^tive immune response; 
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identifying a plurality of biological processes related to the data, the plurality of 
biological processes defining at least one portion of the biological state of the ad£^tive 
unmune response; and 

combining the plurality of biological processes to form an analytical model of the 
adaptive immune response in the context of a peripheral tissue environment and a lymphoid 
tissue environment. 

43. The method of claim 42, wherein at least one biological process ftom the plurality of 
biological processes is associated with recmitment of immune cells into the peripheral 
tissue environment. 

44. The method of claim 43, wherein the inomune cells are blood dendritic cells and blood 
monocytes. 

45. The me&od of claun 44, wherein ttie plurality of biological processes are combined 
' so tibat fhe peripheral tissue enviionmrat is modeled with preferential recruitment of Ae 

biood dendritic cells over the blood monocytes. 

46. The method of claim 42, wherein the biolo^cal state of the adaptive immune 
response is a biological state of an acute response. 

47. The method ofclaim 42, wherein the biological state of flie adaptive immune 
response is a biological state of a chronic response. 

48. The method of claim 42, furOier conq>rising: 

producing an analytical representation of si biological attribute assodated with the 

adaptive immune response; 

•comparing the analytical representation of the biological attribute with a 
corresponding biological attribute in a reference pattern of the adaptive ixzmiune response; 
and 

identifying the analytical model as a valid model of tiie adq>tive immune response if 
the analytical representation of the biological attribute is substantially consistent with the 
biological attribute associated with the reference pattern of fhe adaptive immune response. 

49. The method of claim 42, wherein &e combining the plurality of biological processes 
includes: 
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forming a first mathematical relation among biological variables associated with a 
first biological process firom the plurality of biological processes; and 

forming a second mathematical relation among biological variables associated with 
the first biological process and biological variables associated with a second biological 
process from the pluraUty of biological processes. 

50. The method of claim 42, wherein at least one biological process fi-om the plurality of 
biological processes is associated with a biological variable that is a therapeutic agent. 

51. A method for developmg a computer model of the biological state of an antigen- 
presenting cell, comprising: 

• identifying data relattng to a plurality of physiological regulatory mechanisms of the 
antigen-presenting cell, the data being associated with at least two fix)m the group of antigen 
processing, migration, maturation, and mediator production of the antigen-presenting cell; 

identifying a plurality of biological processes related to the data, the plurality of 
biological processes defining at least one portion of tibie role of the antigen-presenting cell 
in an adaptive immune response; and 

combining the plurality of biological processes to form a simulation of the 
fbnctioning of the antigen-presenting cell in context of the adaptive immune response. 

52. ;The method of claim 51, wherein at least one biological process from the plurality of 
biological processes is associated with a biological variable that is a therapeutic agent 

53. The method of claim 51, wherein the antigen-presenting cell is a dendritic cell. 

54. The method of claim 53, wherein the dendritic cell is a myeloid dendritic cell. 

55. The method of claim 51, wherein at least one biological process fix>m the plurality of 
biological processes is associated with a differential response of lymphocytes to antigen 
based on the maturational state of the antigen-presenting cell. 
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FIG. 4 
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